Zi Gauge Theory of Electron Fractionalization in Strongly Correlated Systems 
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We develop a new theoretical framework for describing and analyzing exotic phases of strongly 
correlated electrons which support excitations with fractional quantum numbers. Starting with 
a class of microscopic models believed to capture much of the essential physics of the cuprate 
superconductors, we derive a new gauge theory - based upon a discrete Ising or Z2 symmetry - which 
interpolates naturally between an antiferromagnetic Mott insulator and a conventional d— wave 
superconductor. We explore the intervening regime, and demonstrate the possible existence of an 
exotic fractionalized insulator - the nodal liquid - as well as various more conventional insulating 
phases exhibiting broken lattice symmetries. A crucial role is played by vortex configurations in the 
Z2 gauge field. Fractionalization is obtained if they are uncondensed. Within the insulating phases, 
the dynamics of these Z2 vortices in two dimensions (2d) is described, after a duality transformation, 
by an Ising model in a transverse field - the Ising spins representing the Z2 vortices. The presence of 
an unusual Berry's phase term in the gauge theory, leads to a doping-dependent "frustration" in the 
dual Ising model, being fully frustrated at half-filling. The Z2 gauge theory is readily generalized 
to a variety of different situations - in particular, it can also describe 3d insulators with fractional 
quantum numbers. We point out that the mechanism of fractionalization for d > 1 is distinct from 
the well-known Id spin-charge separation. Other interesting results include a description of an exotic 
fractionalized superconductor in two or higher dimensions. 



I. INTRODUCTION 

Strongly interacting many-electron systems in low di- 
mensions can exhibit exotic properties, most notably the 
presence of excitations with fractional quantum numbers. 
In these instances the electron is "fractionalized" - ef- 
fectively splintered into consitucnts which essentially be- 
have as free particles. The classic example is the one- 
dimensional (Id) interacting electron gas H, which ex- 
hibits many anomalous properties such as the separa- 
tion of the spin and the charge of the electron. Elec- 
tron "fractionalization" is also predicted to occur in 2d 
systems in very strong magnetic fields that exhibit the 
fractional quantum Hall effect H. Recent experiments 
have given strong supporting evidence of fractionaliza- 
tion both in quantum Hall systems and in carbon nan- 
otubes B. Motivated by these examples, several authors 
have proposed the possibility of electron fractionaliza- 
tion in various other experimental systems. Perhaps the 
most tantalizing is the suggestion by P.W. Anderson g] 
of "spin-charge separation" in cuprate high- Tc materials. 
However, this suggestion is currently surrounded with 
considerable controversy, in part because the Id electron 
gas and the fractional quantum Hall effect appear to be 
rather special situations which do not readily generalize. 
Indeed, in Id the Fermi liquid breaks down even at weak 
coupling and in the quantum Hall regime the kinetic en- 
ergy is strongly quenched by a time reversal breaking 
magnetic field. 

In this paper, we will explore theoretically the possi- 
bility of electron fractionalization in strongly correlated 
systems in spatial dimensions d > 1 in the presence 
of time reversal symmetry. Our primary motivation is 



the cuprates, though we expect our results to be of sig- 
nificance to a variety of other strongly interacting sys- 
tems. Early attempts lq-@] to implement theoretically 
Anderson's suggestion of 2d spin-charge separation typ- 
ically started with either a quantum spin model or the 
t-J model. Slave boson/fermion representations of the 
spin and electron operators were employed to obtain a 
mean field "saddle-point" exhibiting spin-charge sepa- 
ration. The slave boson/fermion representation intro- 
duces a gauge symmetry - f/(l) in the simplest formu- 
lations - and requires inclusion of a corresponding gauge 
field. Fluctuations about the mean field theory lead to 
a strongly interacting gauge theory about which very lit- 
tle is reliably known. It is then quite difficult to reach 
any definitive conclusions about the true low energy be- 
haviour - in particular whether spin-charge separation 
survives beyond the mean field level. An alternate more 
recent approach ^j^ , describes strongly correlated elec- 
tron systems in 2d in a dual language where the vortices 
in the many-electron wavefunction are the fundamental 
degrees of freedom. In this approach, insulating phases 
can be obtained by condensing vortices. Fractionalized 
insulators arise upon condensing pairs of vortices. 

In this work we introduce a new gauge theory approach 
which enables us to reliably address issues of fraction- 
alization. In contrast to the slave boson/fermion rep- 
resentation, our gauge symmetry is discrete -in fact, an 
Ising or Z2 gauge symmetry. This has several advantages. 
Firstly, gauge theories with discrete symmetry are much 
simpler to analyze than those with continuous symme- 
tries nw, so that it is possible for us to make definitive 
statements about low energy physics. But in addition, 
the pure Z2 gauge theory in 2 -I- 1 space-time dimensions 



is dual to the 3d classical Ising model, which implies the 
existence of two distinct quantum phases |12|. In one 
of these two phases "charges" are deconfined, in marked 
contrast to the pure 2 + 1 dimensional U{1) gauge the- 
ory which is always in a confining phase [13[. The pres- 
ence of deconfinement allows us to demonstrate the exis- 
tence of insulating phases exhibiting electron fractional- 
ization, and to describe their basic properties. Remark- 
ably, fractionalization in our Z2 gauge theory approach is 
physically equivalent to vortex pairing in the earlier dual 
formulation P,|l^. We demonstrate this equivalence by 
combining the standard boson-vortex duality ImI with 
the Ising duality mentioned above. 

In addition to the fractionalized phases, our approach 
allows us to readily access the more conventional con- 
fined phases, and the concomitant confinement transi- 
tions. Furthermore, the Z2 gauge theory can be readily 
generalized to describe a variety of different situations - 
arbitrary spatial dimensions, spin-rotation non-invariant 
systems, etc. Some of these generalizations are explored 
towards the end of the paper. For the most part, we con- 
centrate on spin-rotation invariant electronic systems in 
2d. An overview and summary of our main results may 
be found at the end of this introductory section. 

In the context oi frustrated quantvLin spin models. Read 
and Sachdev [|5[ have demonstrated the possibility of 
disordered phases with fractionalization of spin. Specif- 
ically, an Sp{2N) antiferromagnet at large N and the 
related quantum dimer model [n6 17 1 were shown to re- 
duce to a Z2 gauge theory when frustration was present. 
In the deconfined phase of the gauge theory free propa- 
gating spinous (spin 1/2 excitations) would be possible. 
Somewhat similarly, in the slave-fermion representation 
of the conventional Heisenberg magnet which introduces 
an SU{2) gauge invariance, X.G. Wen [Q proposed ob- 
taining fractionalization of spin by pairing and condens- 
ing pairs of spinous. This reduces the gauge symmetry 
down to Z2. In contrast, we show explicitly that the con- 
ventional Heisenberg spin model can be directly written 
as a Z2 gauge theory coupled to fcrmionic spinous, even 
in the absence of any frustration. The key observation is 
that with fcrmionic spinous, the local constraint of sin- 
gle occupancy is equivalent to the constraint of an odd 
number of fermions per site. This latter constraint can 
be implemented with a discrete Z2 gauge field. Such a 
Z2 gauge description may also be obtained with the Ma- 
jorona fermion representation of Heisenberg spins ||l9| . 

The basic physics underlying our description of elec- 
tron fractionalization is perhaps most readily understood 
in d = 2. At the heart of quantum mechanics is wave- 
particle dualism. For a many-body system of interacting 
bosons (with charge Qe, say) this dualism implies that 
in addition to the conventional "particle" framework, a 
description developed in terms of wavefunctions is pos- 
sible. In 2d this dual wave description focusses on point 
like singularities in the phase of the complex wavefunc- 



tion - the familiar vortices with circulation quantized in 
units oi Qy. A fundamental property of such vortices is 
that the product of their quantum of circulation and the 
particle charge is a constant. 



QeQv = he. 



(1) 



It is this simple identity which underlies the two known 
examples of "fractionalization" in two-dimensions, and is 
at the heart of the Z2 gauge theory developed in this pa- 
per. In a ECS superconductor, the pairing of electrons to 
form a Cooper pair with charge Qe — 2e, implies a "halv- 
ing" of the flux quantum, Qy = ^ (hc/e) - which is tanta- 
mount to "vortex fractionalization" . The second example 
of 2d fractionalization occurs in the fractional quantum 
Hall effect g]. In the v = 1/3 state three vortices bind 
to each electron forming a "composite boson" with total 
circulation Qy — 3{hc/e), which then condenses. The 
above identity implies the existence of topological exci- 
tations in this condensate with electrical charge ^e - the 
celebrated Laughlin quasiparticles. 

The route to electron fractionalization that we explore 
in this paper is physieally equivalent to a pairing ofvor- 
tiees, precisely as in earlier work by Balents et. al. P,[lO[. 
But the mathematical implementation is rather different. 
Balents et. al. argued that a pairing and condensation 
of conventional Qy — hc/2e BCS vortices in a singlet 
superconductor, results in an exotic fractionalized insu- 
lator. As Eqn. n^ demonstrates, this insulator should sup- 
port spinless charge e excitations. Our analysis begins by 
noting that such an excitation can be thought of as "one- 
half of a Cooper pair. We implement this fractionaliza- 
tion by formally re-expressing the Cooper pair creation 
operator as the product of two "chargon" operators, 6^, 
each creating a spinless, charge e boson. This change of 
variables introduces a local Z2 symmetry, since it is possi- 
ble to change the sign of 6^ on any given lattice site while 
leaving the Cooper pair operator invariant. This is the 
origin of a local Ising, or Z2, gauge symmetry - described 
mathematically in terms of a Z2 gauge field. In the exotic 
fractionalized insulator, there are strange gapped excita- 
tions which are vortices in the Z2 gauge field. These 
excitations - which we refer to as "visons" because they 
can be represented in terms of Ising spins - are the rem- 
nant of the unpaired hc/2e BCS vortices, which survive 
in the fractionalized insulator. As we shall see, when 
the visons condense they drive "confinement" , thereby 
destroying fractionalization. These visons will play an 
absolutely central role throughout this paper, since any 
insulator with gapped visons is necessarily fractionalized. 

Motivated by the cuprate superconductors, we will fo- 
cus on a particular class of microscopic lattice models 
designed to capture much of the physics believed essen- 
tial to these materials. (Our description of fractional- 
ization is, however, more general and is not restricted 
to these models.) The models describe electrons hop- 



ping on a lattice with inclusion of strong spin and pair- 
ing fluctuations, and are quite similar to models intro- 
duced and analyzed numerically by Assaad et. al. |2^] 
and to models considered more recently by Balents et. 
al. 1^. Many microscopic models of the cuprates, such as 
the t — J model, incorporate spin fluctuations from the 
outset. Our reasons for similarly incorporating "micro- 
scopic" pairing fluctuations are two fold. Firstly, as the 
superconducting phase is a well-established and reason- 
ably well-understood part of the high-Tc phase diagram - 
just like the antiferromagnet - it serves as a useful point 
of departure to access more puzzling regions of the phase 
diagram. This point of view was also advocated in Ref. 
[£l|. But there are also more microscopic reasons to in- 
clude pairing fluctuations from the outset. In particular, 
as emphasized, for instance in Ref. g] , a spin-spin inter- 
action term as in the t — J model can be suggestively 
rewritten in terms of electron operators as. 



iJ'p ' iJt>' 



-^(CrTCr'i - cl^cl,^){h.c) + -prPr'- 



(2) 



term is an attractive pairing interaction in the d^^^ _j,2 (or 
extended-s) wave channel. As in ECS theory, this inter- 
action may be decoupled (in a functional integral) with 
a complex auxiliary pair field rjij as 



E 



J\27]r 
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'{Cr^Cr'l — CrlCj.'^) + C.C.]. (3) 



Here a^-r' = +1 for bonds along the x-direction, and 
equals —1 for bonds along the y-direction. With (?/) ^ 0, 
this corresponds to a superconducting phase with d^2_y2 
symmetry. But more generally, ry can be decomposed 
into an amplitude and a phase, 77 = Ae**^. Ignoring fluc- 
tuations in the amplitude leads to a model of the type 
we consider below, with local fluctuating d— wave pairing 
correlations. 

Further motivation for inclusion of such pairing fluc- 
tuations is provided by resonating valence bond (RVB) 
ideas [p|j22[. The wavefunction for an RVB Mott insu- 
lator can be obtained from the wavefunction of a super- 
conductor by Gutzwiller projecting into a subspace with 
exactly one electron per site. Some mean field theories 
of the RVB state are equivalent to starting out with just 
the superconducting wavefunction. Gauge field fluctua- 
tions about the mean field solution are supposed to carry 
out this highly non-trivial projection, and destroy the su- 
perconductivity. A natural physical route to achieve this 
end is to include strong phase fluctuations of the mean 
fleld order parameter. Indeed, a recent preprint |23| ar- 
gues that fluctuations about the mean field theory of the 
d-wave RVB state pj] are formally equivalent to a theory 
of a phase-fluctuating d-wave superconductor. 

With these motivations, we consider generalized Hub- 
bard type models of the form 



H = Ho + Hj + Hu + Ha 

with , 

Hq = -t^^ cl^Cr'a + h.C, 
{rr') 

Hj ^ J / Sr • Sr' , 

{rr') 

i/„ = ^u(iV,-7Vo)', 

r 

with the local d— wave pair field defined as, 



(4) 

(5) 

(6) 

(7) 
(8) 

(9) 



Here, Cra denotes an electron operator at site r of (say) 
a 2d square lattice with spin polarization a =tii- The 
electron density and spin operators are the usual bilin- 



with pr ~ c].Cr. For antiferromagnetic exchange the first ears: 



Pr 



= rt 



and Sr = kc].crcr with cr a vector 



2 r 



of Pauli matrices. The term i/„ is an on-site repul- 
sion. Strong local tendencies for d^2_y2 pairing are in- 
corporated through the term H^- In the definition of 
Pr in Eqn. BI the summation is over the four nearest 
neighbours of the site r and A^r' — A for bonds along 
the x-direction, and A^r' == —A for bonds along the y- 
direction. With this choice, the operator pr destroys a 
dx2-y2 pair of electrons centered at the site r. 

As discussed above, this anomalous term can be ob- 
tained by decoupling a local spin-exchange interaction - 
which is attractive in the d-wave pairing channel - with a 
complex Hubbard Stratanovich field. Here, we keep the 
amplitude A fixed, but include (quantum) fiuctuations of 
the local pair field phase, Lpr- This phase is canonically 
conjugate to the Cooper pair number operator, Ur'. 



[ipr 



— iSr 



(10) 



Due to the anomalous term in H^, the two densities pr 
and Ur are not separately conserved. The conserved elec- 
trical charge density is simply the sum of the Cooper pair 
and electron densities, 



Nr = 2nr + Pr 



(11) 



It is this total density that enters into the local on-site 
Hubbard interaction term. The c-number A^o plays the 
role of a chemical potential, determining the overall elec- 
trical density. 

This Hamiltonian describes interacting electrons in a 
system with strong local pairing and spin fluctuations. 
Since tpr is a dynamical quantum field, these pairing 
fluctuations do not necessarily lead to a superconduct- 
ing ground state. In addition to the pairing interaction 



terms, the above Hamiltonian includes interactions in the 
spin singlet (u) and spin triplet (J) particle/hole chan- 
nels. The Hamiltonian retains the important global sym- 
metries, corresponding to conservation of spin and elec- 
trical charge. It is worth emphasizing that the theoret- 
ical description of electron fractionalization that we de- 
velop below is not in the least restricted to this particular 
Hamiltonian. 



A. Overview 

Due to the length of this paper, we first provide a brief 
synopsis of our approach and of the key results. We start 
with the observation of Kivelson and Rokhsar [g5[ that, 
in an appropriate sense, a (singlet) superconductor al- 
ready has separation of spin and charge. If one imagines 
inserting an electron into the bulk of a superconductor, 
its charge gets screened out by the condensate to leave 
behind a neutral spin-carrying excitation - a "spinon" . 
A mathematical implementation of this idea [pl| essen- 
tially amounts to binding half of a Cooper pair to an 
electron to produce a neutral spinon. Following these 
ideas, we first split the Cooper pair operator into two 
pieces, each piece creating an excitation with charge e 
but spin zero. These are the same quantum numbers as 
the "holon" . But since this object seems to be defined 
rather differently, and in any event is not directly tied to 
the doping of a Mott insulator, we prefer to refer to it 
as a "chargon" . The square of the chargon operator cre- 
ates the Cooper pair. Next, we define a neutral "spinon" 
operator by multiplying together the chargon and elec- 
tron operators. Changing variables from the electrons 
and Cooper pairs to chargons and spinous introduces a 
degree of redundancy in the description. Specifically, all 
physical observables are invariant under a local change 
in the sign of the spinon and chargon operators. This 
implies that the resulting theory must have a local Z2 
gauge invariance. 

In Section y, we carefully re-express the above model 
in terms of the chargon and spinon operators, paying spe- 
cial attention to the local Z2 gauge symmetry. Following 
techniques familiar from slave boson/fermion theories, we 
derive an action in terms of the chargon and spinon fields 
coupled to a fluctuating Z2 gauge field. This takes the 
form 



time lattice labelled by i, j, The spin is carried by the 

(Grassmann-valued) spinon fields, fi and /i, also living 
on the lattice sites. The chargon and spinon fields are 
"minimally coupled" to an Ising Z2 gauge field CTy = ±1 
living on the links of the space-time lattice. The form of 
the charge and spin actions, Sc and Ss, could have been 
guessed on symmetry grounds (the global charge t/(l), 
the global spin SU(2) and the local Z2 gauge symmetry), 
but the derivation in Section || shows the presence of an 
additional term Sb- This is a "Berry phase" term that 
takes the form 



Se^-i Yl Noi'^T^hj-^il-a.,^ 



(16) 



Here t refers to the time direction, and hj is an integer 
on each temporal link defined in terms of the (p and a 
fields as, 



Int 



2tt 



(17) 



with $jj the gauge invariant phase difference across the 
temporal link: 



$,; 



^i+ 2(l~'^y)- 



(18) 



The symbol Int refers to the integer part. The Berry 
phase term simplifies considerably for integer Nq. For 
even integer Nq, we simply have e~^ — 1, while for odd 
integer A^o, 



,~Sb _ 



n 



No odd. 



(19) 



■,,J=l-T 



A rough estimate of the dimensionless couplings 
tcV^ , t^ in terms of the parameters i, u, J, A of the orig- 
inal microscopic Hamiltonian may be obtained in the 
physically interesting limit of large u and small t near 
half-filling: 



J 



tc\ f- 



A 
—i 
t 



(20) 



S ^ Sc + Ss + Sb , 
Sc = -tc ^ (Jij {b*bj + c.c), 



(12) 
(13) 



(ij) 



Ss^ -Y^ (Jij {ttjfiafja + ttjli^lji + C-C) - Y fio^fio.- (l^) 



(u) 



We will, however, regard these coupling strengths as phe- 
nomenological input parameters for the Z2 gauge theory. 
A great deal of physics is contained in the simple- 
looking action Eqn. n3. Consider varying the dimen- 
sionless chargon coupling, ic, which represents the degree 
of charge fluctuations, and for simplicity specializing to 
half- filling with Nq = 1. Surprisingly, in the limit of 
vanishing chargon coupling, t^ — 0, the full Z2 gauge 
theory action can be shown to be formally equivalent 



(15) 



Here Sc describes the charge dynamics with hi = e "^' 
the chargon field defined on a d -|- 1 dimensional space- 



(See Section IV) to the Heisenberg antiferromagnetic spin 
model. Increasing tc from zero introduces charge fluctu- 
ations into the Heisenberg model. In the limit of large 
tc, the chargons will condense resulting in a conventional 



da:2-y2 superconductor. Thus, the above Z2 gauge theory 
action has the remarkable property of interpolating be- 
tween the Heisenberg antiferromagnet in one limit and a 
dx2_y2 superconductor in the opposite limit. Determin- 
ing the properties of this model in the intervening regime 
(with tc of order one) is an extremely interesting question 
in the context of the cuprate materials, and will be one 
of the prime focuses of our analysis. Specifically, within 
the present Z2 gauge theory we will explore the different 
possible routes between these two limits (which depend 
on the parameters in the action). Most importantly, for 
certain parameter regimes we will demonstrate the pos- 
sibility of obtaining an exotic fractionalized insulating 
phase - dubbed the nodal liquid in previous work |2l|] - 
intervening between the antiferromagnet and the (i^2_j,2 
superconductor. For other parameter regimes, a number 
of conventional insulating phases (z.e, with no fraction- 
alization) are accessible, including various phases with 
spin-Peierls and/or charge order. 

To gain a simple understanding of these results it is ex- 
tremely convenient to integrate out the chargons to give 
an effective action depending only on the spinous and the 
Z2 gauge field a. This is legitimate provided the char- 
gons are gapped, as they will be in all of the insulating 
phases (with iVo = 1). The most important effect of this 
integration will be to generate a "kinetic" term for the 
Z2 gauge field a: 



S^ = -KJ2[l['^^. 



(21) 



Here, the product is of the Z2 gauge fields around an 
elementary plaquette of the space-time lattice, and this 
product is then summed over all plaquettes. Clearly, Sa 
is the direct Ising analog of the F^^ term which enters 
the Lagrangian of ordinary C/(l) electromagnetism. The 
value of the parameter K is determined by the chargon 
coupling, increasing monotonically with tc- The full effec- 
tive action appropriate to the insulating phases is simply. 



S — Ss + Sa -\- Si 



(22) 



Since the onset of superconductivity will occur at some 
critical value of order one, t* « 1, the validity of the ef- 
fective action requires tc < t*. Near this limit, but on 
the insulating side, K will also be of order one. 

There are several limits in which the properties of this 
effective action may be reliably analysed. A schematic 
phase diagram is shown in Fig. n. As mentioned above, 
with K = tc = the action describes the Heisenberg 
antiferromagnet, which in 2d exhibits Neel long-ranged 
order at zero temperature. The opposite limit of large K 
is far more interesting, though. Indeed, when K = 00, 
fluctuations of the Z2 gauge field (Ty are frozen, and one 
can set (Ty ~ 1 on all the links. This results in a phase 
with deconfined spinous propagating freely with the gap- 
less "d-wave" dispersion - the "nodal liquid" . Similarly, 



the chargons are also deconfined, existing as gapped exci- 
tations in this insulating phase. The nodal liquid is thus 
a genuinely "fractionalized" insulator, within which the 
electron has splintered into two pieces that propagate in- 
dependently. On reducing K from 00, the nodal liquid 
continues to be stable until a certain critical value Kc 
of order one, where the gauge field undergoes a confine- 
ment transition. For K < Kc the chargons and spinous 
are no longer legitimate excitations, but rather are con- 
fined together to form the electron (or other composites 
built from the electron such as magnons or Cooper pairs) . 
This corresponds to a conventional insulating phase. As 
we argue in Section |V|, the confinement transition is ac- 
companied by a breaking of translational symmetry lead- 
ing to spin-Peierls order - at least for small spinon cou- 
plings V^ ,t^ . This may be understood from the limit 



when t^ ,t = 0. Then, as we show in Section IV, we are 



left with a pure Z2 gauge theory with the Berry phase 
term Sb which is exactly dual to the fully frustrated Ising 
model in a transverse magnetic field. Ordering the Ising 
spins in this dual global Ising model leads to confine- 
ment. Physically, the Ising spins represent vortices in 
the Z2 gauge field - namely, the vison excitations men- 
tioned in the previous subsection. This same model also 
arose in the studies of Sachdev and coworkers [|l6|,0 on 
frustrated large N quantum antiferromagnets. Numeri- 
cal studies [|l6| show that the ordering in the Ising model 
is accompanied by breaking of translational symmetry. 
The nature of the confined phase(s) at large spinon cou- 
pling remains uncertain at present. 




FIG. 1. Schematic zero temperature phase diagram of the 
insulating phases showing the three limits mentioned in the 
text. The horizontal axis measures the strength of the cou- 
pling K obtained by integrating out the chargons. The ver- 
tical axis is a measure of the spinon couplings 1^,1^. Here 
AF denotes the Heisenberg antiferromagnet, SP denotes an 
insulator with broken translational and rotational invariances 
such as a spin-Peierls state, and NL denotes the nodal liquid 
with fractionalized excitations 



These results demonstrate the possibility of two alter- 
nate routes between an antiferromagnet and a d— wave 
superconductor. In one instance, as the chargon hopping 
tc is increased towards the critical value for the onset 
of superconductivity t*, the value of the parameter K 
stays smaller than the critical value for deconfinement, 
Kc- In this case, all of the insulating phases preceding 
the superconductor are "conventional" , with confinement 
of chargons and spinous. Alternately, if K exceeds Kc be- 
fore the transition into the superconductor, the fraction- 
alized nodal liquid phase will occur - sandwiched between 
the d— wave superconductor and a conventional insulator. 
Since both the superconducting and the deconfinement 
transitions occur when tc (and hence K) is of order one, 
the deconfinement boundary is expected to be "near" the 
onset of superconductity It is thus difficult to ascertain 
which of these two scenarios will be realized. The precise 
phase diagram interpolating between the antiferomagent 
and superconductor will likely depend sensitively on var- 
ious microscopic details. 

Considerable further insight is provided into the mech- 
anism of electron fractionalization in an alternate dual 
formulation in which we trade the chargon fields for the 
hc/2e vortices which occur in a conventional supercon- 
ductor. In Appendix H, we show how this may be done 
following standard duality transformations for the classi- 
cal three dimensional XY model. Starting with the full 
Z2 gauge theory in Eqn. |lj, the resulting dual theory 
is a lattice action for the hc/2e vortices coupled to the 
spinous. The vortices see a fluctuating f/(l) gauge field 
a whose circulation is the total electrical three current. 
Further, the hc/2e vortices have a long-ranged statistical 
interaction with the spinous: When a spinon encircles 
such a vortex, its wavefunction acquires a phase of tt. In 
the present formulation, a tt flux of the Z2 gauge fleld a 
is effectively attached to each vortex. As the spinous are 
minimally coupled to u, they acquire the expected phase 
of TT upon encircling each vortex. Mathematically, this 
flux attachment is implemented by an analog of a Chern- 
Simons term for the Ising group. Quite remarkably, this 
Ising Chern-Simons term emerges automatically from the 
duality transformation in Appendix ^ 

This dual representation of the Z2 gauge theory is in 
fact essentially identical to the vortex field theory intro- 
duced in Ref. [g| on a phenomenological basis starting 
with a BCS superconductor. In that work, the statisti- 
cal interaction between spinous and vortices was put in 
by hand, employing a U{\) Chern-Simons terms to at- 
tach flux to the spin of the spinous. An advantage of 
the Ising Chern-Simons terms is that it does not break 
spin-rotational invariance, and in fact is possible even for 
spinless electrons. Moreover, it enables the description of 
an exotic superconducing phase in which the Ising flux 
de-attaches from the vortices (see below). In this dual 
description, superconducting phases correspond to vor- 
tex vacuua, while insulating phases correspond to vortex 



condensates. Simply condensing the hc/2e vortices leads 
to conflned insulating phases. Accessing deconfined in- 
sulating phases requires condensation of paired vortices, 
without condensation of single ones. In this way one ob- 
tains an alternate dual description of the fractionalized 
nodal liquid. The Z2 gauge theory formulation suggests 
a mechanism for such vortex pairing: Since the chargons 
also have a long-ranged statistical interaction with hc/2e 
vortices their motion is "frustrated" in the presence of 
such vortices. Pairing the vortices reduces this frustra- 
tion, allowing the charge to propagate more easily, and 
lowering the kinetic energy. 

Superconducting phases are readily accessed in either 
the Z2 gauge theory "particle" formulation of Eqn. |l^, 
or its dual vortex counterpart. In the particle formula- 
tion, when tc becomes large and the charge e chargons 
condense, the result is a (1^2 _yi superconductor - denoted 
dSC. This superconductor is conventional, perhaps sur- 
prising since BCS theory involves the condensate of a 
charge 2e Cooper pair. But as we demonstrate in Sec- 
tion [I], the chargon condensate supports hc/2e vortices. 



and shares all other properties with a conventional BCS 
superconductor. It is interesting to ask if it is possible 
to have a superconductor where the chargon pairs have 
condensed, while the single chargons have not. Such a 
superconductor, which we label dSC* , can be readily de- 
scribed with the present Z2 gauge theory formulation. 
As detailed in Sec. VII, dSC* is a truly exotic supercon- 



ducting phase with many unusual properties. 

The Z2 gauge theory is readily generalized to a wide 
variety of other situations. In particular, the "particle" 



formulation of Eqn. 12 is valid in any spatial dimension. 



In 3d there again exist fractionalized insulating phases 
(and of course confined ones) which can be a ccessed by 
the theory. Remarkably, as we argue in Section VIII E , in 
contrast to the 2d case, a fractionalized insulator in 3d 
exists as a distinct finite temperature phase, separated 
by a classical phase transition from the high tempera- 
ture limit. For an anisotropic layered three dimensional 
material, it is also possible to have another 3d fraction- 
alized phase consisting of weakly coupled 2d phases, but 
this phase is destroyed by thermal fluctuations. It is also 
noteworthy that the Z2 gauge theory formulation seems 
incapable of describing fractionalization in Id. This indi- 
cates that the "solitonic" mechanism of fractionalization 
in d = 1, is qualitatively different than "vortex pairing" 
which describes fractionalization in higher dimensions. 

We conclude this section with an outline of the rest of 
the paper. Section O contains the formal derivation of 
the Z2 gauge theory from the microscopic models. For 
ease of presentation, and as it is simpler, we will first 
provide the technical details of the derivation for situa- 
tions with local s-wave pairing. In Appendix pi we show 
how situations with d^2_y2 pairing, the case of interest 
for the cuprates, can be readily handled. We next de- 
scribe, in Section fn, the physics of fractionalization and 



confinement in the simplest possible context - that of s- 
wave pairing with an even number of electrons per unit 
cell. We then consider, in Section |^, the more interest- 
ing situation of d-wave pairing with an odd number of 
electrons per unit cell. Section |^ formulates and devel- 
ops the dual description in terms of vortices. The results 
of Section IV are reobtaincd in this representation. We 



then move on in Section VI to show how doping away 



from half-filling may be incorporated into the formalism. 



In Section VII , we discuss the possibility of other exotic 



fractionalized phases, in particular the superconductor 
SC* mentioned above, in both the particle and vortex 
formulations. Section VII] discusses various generaliza- 
tions of the theory, including spatial dimensions other 
than two, finite temperature, and situations with no spin 
rotational invariance. We also briefly discuss a useful 
analogy with Z2 lattice gauge theories of classical ne- 
matic systems. In Section IX, we discuss the relationship 
between this work, and several other previous approaches 
to fractionalization in strongly correlated systems. Con- 
tact will be made, when possible, with the earlier dual 
vortex descriptions of the nodal liquid, and with the slave 
boson/fermion approaches. Section M contains a discus- 
sion of the experimental signatures of the various novel 
phases described in earlier sections. We conclude with a 
summary of our main results. Various appendices con- 
tain technical details not presented in the main text. 



II. MODELS AND Z2 GAUGE THEORY 

To describe our techniques in the simplest possible con- 
text, we will start with a microscopic model that has local 
s-wave pairing correlations. This can be readily general- 
ized to other symmetries such as d-wave (see the end of 
this Section, and Appendix ^). Of course, with strong 
local on-site repulsion (positive u above) d-wave pair- 
ing fluctuations are presumably more energetically vi- 
able, and also of central interest in the context of cuprate 
superconductivity. 

Consider then a generalized Hubbard type model: 



H^Ho+Hu + Hj + Ha, 



with 



{rr') 



t X! 4aCr'a + h.C, 



No)', 



Hj = j'^ [Sr ■ Sr' + -prPr'], 

(rr'} 

TJa = A ^(e^'^'-c^-rc^i -t- h.c). 



(23) 

(24) 
(25) 
(26) 
(27) 



As earlier, Cra denotes an electron operator at site r with 
spin a and the electron density and spin operators are 



the usual bihnears: pr = c\.a^ra and Sr = iclcrcr- This 
Hamiltonian is essentially the same as Eqn. ^ in the intro- 
duction, except that it has local s— wave pairing rather 
than d— wave, and we have added a term proportional 
to prPr' in Hj. These modifications have been made to 
both simplify the derivation and the subsequent analy- 
sis of the Z2 gauge theory. We return later to the more 
physically interesting case of local d— wave pairing. 

Here, (pr is the phase of a local s-wave Cooper pair field 
and is canonically conjugate to the Cooper pair number 
operator, n^: [(/9r,nr'] = iSrr'- As before, since ipr is a 
dynamical quantum field, these pairing fluctuations do 
not necessarily lead to a superconducting ground state. 
The conserved electrical charge density is the sum of the 
Cooper pair and electron densities 



Nr = 2nr + Pr- 



A. Split the Cooper pair 



(28) 



We now proceed to split the Cooper pair into two 
pieces. Consider an operator br defined as, 



bl 






iipr/2 



J4>T 



(29) 



with Sr = ±1 an Ising "spin" variable. With this defini- 
tion the new field, 



2 2^ 



(30) 



can be treated as a phase - lying in the interval zero to 27r, 
with br invariant under the transformation: ipr -^ ipr+2T: 
and Sr — > —Sr- The square of bl creates a Cooper pair. 



p«v 



-- - (bir 



(31) 



so that bl creates a spinless excitation with charge e, 
essentially one-half of a Cooper pair. We refer to this 
operator as a "chargon" operator. 

In order to separate out the charge and spin degrees 
of freedom it will be extremely useful to define an elec- 
trically neutral but spin carrying fcrmion operator (a 
"spinon"): 



f 1^ — /) r^ 

Jra ~ "r<~-rc 



(32) 



This operator carries the spin of the electron, but is elec- 
trically neutral as verified by noting that it commutes 
with the total electrical charge density Nr- On the other 
hand, the chargon is electrically charged, and its phase is 
canonically conjugate to the total electrical charge den- 
sity 



[(j)r,Nr 



iSr 



(33) 



At this stage it is legitimate to implement an operator 
change of variables in the full Hamiltonian, replacing the 
electron and Cooper pair operators {(p,n,c,c') by char- 
gons and spinons {(f), N, f,f^). This gives, 



with 



H = Ho + Hu + Hj + Ha, 



Ho = -t ^ blK'flafr'a + h.C, 

{rr') 



HA = Aj2iMfri + h.c.), 



(34) 



(35) 



(36) 



with Hu unchanged and Hj of the same form as in Eqn. 
Pq but with spinon operators replacing the electron op- 
erators: Pr = flafra and Sr = \!}.a-fr- 

There are several extremely important points to stress 
about this seemingly inoccuous change of variables. 
Firstly, one can change the sign of both the chargon and 
spinon operators on any given site r. 



fr 



Jro 



(37) 



without affecting the original Cooper pair or electron 
operators. This implies that quite generally the trans- 
formed Hamiltonian must also be invariant under this 
local Ising Z2 symmetry - as can be readily checked in 
Eqns. p^ and pq. As we shall shortly see, in a path inte- 
gral formulation this local Z2 symmetry will be manifest 
in terms of an Z2 gauge field. Secondly, because of this 
redundancy introduced in the change of variables, a con- 
straint must be imposed on the Hilbert space spanned by 
the spinon and chargon operators. 

To understand the origin of this constraint, consider 
first the Hilbert space of the original Hamiltonian. In 
a number-diagonal basis, the Hilbert space on each site 
r is a direct product of states with an arbitrary integer 
number of Cooper pairs (n^) and the four electron states 
consistent with Pauli - empty, doubly occupied or singly 
occupied with an electron of either spin. Since the char- 
gon has only one-half the charge of the Cooper pair, the 
full Hilbert space spanned by the chargon and spinon 
operators is actually twice as large, and it is essential to 
project down into the physical Hilbert space of electrons 
and Cooper pairs. From Eqn. ^, it is clear that this can 
be achieved by imposing a constraint that the sum (or 
difference) of the number of chargons (Nr) and spinons 
{Pr = frafra) On cach sitc is an even integer: 



(-1) 



N^+P,^ _ 



= 1. 



(38) 



This implies, for example, that a site with a single char- 
gon but no spinon is unphysical and forbidden, whereas 
a spinon and chargon together (an electron) is allowed. 



B. Path Integral and Z2 gauge Theory 

The most convenient way to implement the constraint 
on the spinon and chargon Hilbert space is in a (Euclid- 
ian) path integral representation of the partition func- 
tion. To this end we define a projection operator. 



r = l[rr 



(39) 



with 



7'. = -[l + (-l)^'-+'"- 






„if(l-<T,)(Ar,,-Hp,,) 



(40) 



which projects into the physical Hilbert space. Here, 
CTr = ±1 is an Ising-like field and pr = frafra- As can be 
verified directly from Eqn. |3J, this projection operator 
commutes with the chargon-spinon Hamiltonian, 



[V,H] = 0, 



(41) 



so that the Hamiltonian does not cause transitions out 
of the physical Hilbert space. 

The partition function can be written as, 



Z = Tr[e-^"r], 



(42) 



where the trace is over the full Hilbert space spanned by 
the chargon and spinon operators (0, N, f, f^). A Euclid- 
ian path integral representation can be obtained as usual 
by splitting the exponential, 



Tr[{e-'"Vf% 



(43) 



with M "time slices" and e = (3/M . Here, we have in- 
serted projection operators into each time slice. Working 
with fermion coherent states and eigenstates of the char- 
gon phase (f), a path integral representation can be readily 
derived - as detailed in Appendix A - giving, 



Z = 



/oo 
Wdjiadfiad(t>r ^ ^ 



(44) 



Ni = -OQ iJi=±l 



where the integration is over Grassman numbers / and / 
and a c-number phase (j) in the interval zero to 27r. Here, 
i — (r, r) runs over the 2 -|- 1-dimensional space time lat- 
tice with r = 1, 2, ..., M time slices. The Euclidian action 
takes the form, 

M 
S^Sl+St + eY,H{Nr,ct)r,Jrfr), (45) 



T = l 



ith 



M 



5/= ^[/.(a.+i/.+i-/.)], (46) 

r,T— 1 

M 

5^ = -* ^ iV.[(/.. -</)._! + |(1- a,)]. (47) 



where 6* = e*'^''. The full Euclidian action then takes 
the form, S = S( + Sf + Sr, with the spatial interactions 
given by, 



Sr^eJ2iHu + HA) + S^ 



(54) 



Here, we have suppressed the explicit r and a subscripts 
on the fields, displaying only the time-slice dependencies. 
As usual, the bosonic phase field and the Ising field both 
have the expected periodic boundary conditions, whereas 
the fermions are anti-periodic: 

0r=M+l = 0r = i; CTm+I = CTi ; /m+1 = -/l- (48) 

Notice that the Ising variables live on the links con- 
necting adjacent time slices, and can thus be correctly 
interpreted as a gauge field. In fact, the Ising field a 
is minimally coupled to both spinous and chargons, as 
the time component of a gauge field. Moreover, the local 
Z2 symmetry of the Hamiltonian in Eqn. |3^, is mani- 
fest in the path integral as a full fledged Ising Z2 gauge 
symmetry: 



. + -(l-eO, (49) 



Jia ^ ^ijia-! Jia ^ ^ijia-i (pi 

together with a transformation of the gauge field, 



(50) 



Here, e^ = ±1, and aij lives on the link connecting two 
"nearest neighbor" space-time lattice points, differing by 
one time slice. 

Our final goal is to beat the model into a form which 
also includes Z2 gauge fields on the spatial links, so that 
space and time end up on more equal footing. Our ap- 
proach follows closely the standard methods pq] em- 
ployed in slave fermion or slave boson treatments of 
Heisenberg magnets. First, we perform a Hubbard- 
Stratanovich decoupling of the spin interaction terms in 
the Euclidian action: 



e-ff. 



n n^^'''''^'^)^^*'-'^'^)^' 



Shs 



(51) 



{rr') 



Shs = e ^ ^[2J|Xrr'P - {JXrr'frafr'a + C.C.)]. (52) 



{rr') 



Here, Xrr' (t) are a set of complex fields which live on 
each of the nearest neighbor spatial links. Next, a simple 
change of variables can be performed which eliminates 
the remaining quartic spinon-chargon interaction, in Hq 
in Eqn. p3: 



Xr 



Xr 



Kbr': 



(53) 



with 

S'x = e ^ 2J\Xrr'f - [Xrr'{'2tb*,br + Jfrafr'a) + C.C.] 



{rr'} 



(55) 



The terms in S^^^ correspond to the hopping of spinous 
and chargons in the presence of a common fluctuating 
gauge field, %, on the near neighbor links. 

Up to this stage, all of the formal manipulations that 
we have performed have been exact, so that the full Eu- 
clidian action gives a faithful representation of the orig- 
inal "microscopic" electron Hamiltonian. But now, fol- 
lowing standard slave fermion/boson techniques, we per- 
form an approximation, treating the functional integral 
over the Hubbard-Stratanovich field, x-i within a sad- 
dle point approximation. [While it might be possible to 
find an appropriate "large-N" generalization of the model 
for which this approximation becomes exact, we do not 
pursue this tack here.] The simplest saddle-point corre- 
sponds to setting all of the link fields equal to a single 
real constant: Xrr' — Xo- The saddle-point value for 
Xo can (in principle) be obtained by integrating out the 
spinous (which are Gaussian) and the chargons (which 
are not). This saddle-point respects two important dis- 
crete symmetries of the model - translational and time- 
reversal invariance. But the saddle-point does noi respect 
the Z2 gauge symmetry in Eqns. p^ and pOJ. This serious 
fiaw can be easily remedied, though, by retaining a par- 
ticular set of fluctuations about the saddle point. The 
simplest choice consistent with the Z2 gauge symmetry 
corresponds to allowing the sign of Xrr' to change, keep- 
ing the magnitude fixed, putting 



Xr 



'Xo- 



(56) 



Here, arr'{T) — ±1 are a set of Ising fields living on 
the spatial links of the space-time lattice. Within this 
restricted manifold the theory consists of chargons and 
spinous hopping on a space-time lattice, minimally cou- 
pled to an Z2 gauge field. 

Hereafter we work under this fixed-magnitude approx- 
imation. Within this approximation the full partition 
function can be expressed as a functional integral. 



/oo 
i[df,^dho.d(t), ^ n E ^"^' 
,■ -1, AT -^ / ■ ■ \ ^ . . LI 



(57) 



Ni = -00 {ij) (Tij=±l 



with Z2 gauge fields Uij living on the near neighbor links 
of the space-time lattice, and 



s — si + Sf + Sq + Su + Sa, 
with 

Jt — 2-^ [jia{o'ijJja — Ji)\, 



(58) with the "Berry phase" term Sb given by 



—i 



St^-l Y. m'l^^-^J+^{^-<^^J)], 



Su^€uJ2{N^-No)^, 

i 

SA=eAY,{MM + kM), 

i 

5o = -e ^ aij{tob*bj + Jofiafja +c.c.), 



(59) 
(60) 
(61) 
(62) 
(63) 



i,j^i+x 



where we have defined tg = 2txo and Jq ~ Jxo- 

Notice that the full action is local in the integers Ni , so 
the summation can be performed independently at each 
space-time point. A straightforward Poisson resumma- 
tion gives 



Y^e~i'-^'i)=eM Y: n$.)], 



(64) 



Ni 



■,,J=l-T 



where $ij = (pi — (pj + ^{\ — Uij) is the gauge invariant 
phase difference along a temporal link. Here, the periodic 
potential V($) is given by 



:,vm _ 



/— — OO 



,- ^[<i--2TTl\^ +iNo(2T,l-'i>) 



(65) 



and we have dropped an overall multiplicative constant. 
In the limit of small eu, the sum over I will be dominated 
by precisely one term which minimizes |$ — 27r/|. This 
occurs for integer I satisfying |$ — 2t:1\ < -k, or equva- 
lently, 



$ 1 

Moreover, for small eu we may approximate 



e 4 



ij(*-2;T0^ 






(66) 



(67) 
(68) 



Within this approximation the sum over I becomes sim- 
ply, 



,y($) 



+ ^cos(-i>)+iNa(2TTl-<^) 



(69) 



with I given by Eqn. pq. We have again dropped an 
overall multiplicative constant. 

The full N sum in the action then leads to 



^e-(^"+^*)=eS'.^-— ^^"•- 



jCOs(ct>i—<pj) — SB 



(70) 



Sb = -iNo J2 (^T^^y-^y)' 
= -iNo Y [27rZy -|(l-fJi 



(71) 
(72) 



*j=*-'r 



In obtaining the last line, we have re-expressed $ij in 
terms of (p and a, and used the /?— periodic boundary 
conditions on p to drop the term involving p>i — (pj . The 
"Berry phase" term is the only term in the action which 
depends on the (average) occupation number per unit 
cell, iVo- It simplifies considerably for integer A^o- For 
even integer Nq, we simply have e~^ = 1, while for odd 
integer Nq, 



n 



iVo odd. 



(73) 



As we shall see, the Berry's phase term will lead to subtle 
yet important differences between Mott insulators with 
odd integer A'o and band insulators with even A^o- 

The Euclidian path integral is only identical to the 
Hamiltonian formulation in the strict e — > limit. But 
since the original lattice Hamiltonian is already an effec- 
tive low energy theory, the time continuum limit which 
involves arbitrarily high energies is not actually of inter- 
est. For these reasons, hereafter we keep e finite, viewing 
it as an inverse "high energy" cutoff in the theory. Since 
the kinetic (i) and interaction (u) energy scales are the 
largest in the theory, it is convenient to choose the value 
of e so that the charge sector of the theory is isotropic on 
the 2 + 1-dimensional space-time lattice. To this end, we 
require that the spatial chargon hopping strength equals 
the temporal one: ^— = 2eto, which implies 



2\/tc\u. 



(74) 



With this choice of e the full Euclidian action reduces 
to a much simpler and more compact form: 



(75) 



(76) 



S — Sc + Ss + Sb 
with 
Sc = ~tc Y, ^viK^j + h.c), 

Ss = Y. -^h'^-ohfj + c.c.) + 5.,(i^/.T/^i + c.c. - /,/,), (77) 

and Sb as defined above. Here, the dimensionless char- 
gon coupling strength is given in terms of the microscopic 
parameters t, u and xo to be 



Ni 



tc — ein — 



(78) 
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The dimensionless spinon coupling along the nearest 
neighbor spatial links is 



if, = e Jn = J 




(79) 



whereas tf. = — 1 along the neighboring temporal links. 
Similarly, the coupling constant for the spinon pairing is 



t^ 



A 



\/8tXoU 



(80) 



As will be shown in Section IV , for the physically interest- 
ing case of d-wave pairing near half-filling, the parameter 
Xo niay be roughly estimated to be 



Xo 



J2 



(81) 



This can be used to obtain rough estimates of the three 
dimensionless coupling constants, tc,t^ and t^. For the 
most part, however, we will treat these couplings as phe- 
nomenological parameters. 

The partition function involves an integration over the 
on-site chargon phase {(pi) and spinon Grassman fields 
ifi: fi)i as well as a summation over the Z2 gauge fields 
{aij = ±1) which live on the nearest neighbor links of the 
Euclidian space time lattice. This "final" form for the 
theory is exceedingly simple, consisting of chargons and 
spinous hopping around, minimally coupled to a dynam- 
ical Z2 gauge field. This form could have essentially been 
guessed just using a knowledge of the field content (char- 
gons and spinous) and the required symmetries; U{1) 
charge conservation, SU{2) spin conservation and the lo- 
cal Z2 gauge symmetry. Perhaps the only subtlety is the 
presence of the term Sb in the action when the filling 
factor A^o is not an even integer. Among the additional 
terms which are allowed by these symmetries, is a field 
strength term for the Z2 gauge field: 



Sa = -Kj^iU'''.^^- 



(82) 



Here, the product denotes the gauge invariant product 
of the Ising fields around an elementary plaquette. This 
Ising field strength is then summed over all space-time 
plaquettes. Clearly, Sa is the direct Ising analog of 
the F^^ term which enters the Lagrangian of ordinary 
C/(l) electromagnetism. Even though not present in the 
derivation presented here, this field strength term will 
be generated upon integrating out the chargon or spinon 
matter fields, as discussed below. 

In Appendix M we show how the above analysis can 
be generalized to the case in which local d-wave pair- 
ing correlations are incorporated from the outset as in 
the Hamiltonian Eqn. H, rather than s-wave as assumed 
above. The derivation of the effective Z2 gauge theory 
proceeds in much the same fashion, and one arrives at the 



same model except with the spinon action given instead 
by 

'-'s — / ^ "'ij V'ijjiajja + ^ijJiT Jjl ' C.C) / ^ Jiajia- 



iv) 



(83) 



Here, t~^ denotes a d-wave pairing amplitude living on 
the nearest neighbor spatial bonds, with amplitude -t-i"^ 
on the X-axis bonds and —f^ along the y-axis bonds. No- 
tice that the Z2 gauge field aij enters here, because the 
d-wave pair field lives on the links. This form exhibits the 
required Ising Z2 gauge symmetry, being invariant under 
the transformation in Eqn. ^. As shown in Section IV, 
a rough estimate of the various coupling constants in this 
case is 



tc 



J 



u 



f'^^ 



f- 



jtc- (84) 



Here i'* and f^ refer only to the spatial couplings. But 
we will once again regard these as phenomenological pa- 
rameters. 



III. FRACTIONALIZATION AND 

CONFINEMENT 

In this section we will analyze some of the phases which 
are described by the Z2 gauge theory model derived in 
Section II. While the Z2 gauge formulation is valid in gen- 
eral dimension, for concreteness and simplicity we spe- 
cialize to two dimensions, generalizing briefly to other 
dimensions in Section VIII A. Moreover, for illustrative 



purposes we focus first on the simplest case with an even 
number of electrons per site (unit cell), and presume the 
presence of local s-wave pairing correlations. As we shall 
see, in this case the model can exhibit a conventional 
band insulator. In Section IV we will turn to the more 



physically interesting situation with an odd number of 
electrons per site. At that stage we will focus on local d- 
wave pairing correlations, which are more tenable in the 
presence of a large positive on-site Hubbard u as well as 
being of direct relevance to the cuprates. Doping away 
from half-filling will be discussed in Section [Vl| . 

With even integer A^o and local s-wave pairing cor- 
relations the full action consists of two contributions, 
S — Sc + Ss, corresponding to the charge and spin sec- 
tors, respectively: 



Sc = -tc ^ cTij {b*bj + c.c), 



(ij) 



Ss^-Y. ^3^n (Mj + c.c.) - J2 f^f^ 



(ij) 



<^E(/'T/a + c.c.). 



(85) 



(86) 



(87) 
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The first term, which describes the dynamics of the char- 
gons, b* = e**^, minimally coupled to an Z2 gauge field, 
exhibits the global U{1) charge conservation symmetry. 
The spinous also carry the Z2 Ising "charge" . Due to the 
s- wave form of the anomalous "pairing" term the spinous, 
which are paired into singlets, should be gapped out. 



1. Correlated "Band" Insulators 

We first consider electrically insulating states. When 
the dimensionless chargon coupling tc is much smaller 
than unity, the chargons cannot propagate at low ener- 
gies and a charge gap results. In this case, with both 
spinous and chargons gapped out, it is possible to inte- 
grate them out from the theory, leaving the Z2 gauge 
field (7 as the only remaining field. This integration will 
generate additional terms in the Lagrangian, depending 
on (J, which will be local in space-time and must also be 
gauge invariant. The most important such term P6[ is 
simply. 



S. = -KY,[Y[a, 



(88) 



which describes a pure Z2 gauge theory. 

Remarkably, this simple gauge theory exhibits a phase 
transition as the coupling K is varied. Indeed, as shown 
originally by Wegner [|l^,0 , the pure Z2 gauge theory in 
3D is dual to the familiar three-dimensional Ising model: 



Sdual = -Kd ^ ViV 
(iJ) 



J' 



(89) 



with Ising spins, Vi = ±1, living on the sites of the dual 
lattice. The dimensionless Ising model coupling, Kd, 
is simply related to K: tanh{Kd) = e~^^ . This form 
shows that the high and low "temperature" phases are 
exchanged under the duality transformation. The details 
of this duality transformation are given in Appendix |C[ 
As emphasized originally by Wilson |2^ , a direct char- 
acterization of the two phases of the pure gauge theory 
is given in terms of the correlator. 



5c = {ll^ 



(90) 



where the average is for the pure gauge theory and the 
product is taken around a closed loop in space-time, de- 
noted C. For K < Kc the Wilson-loop satisfies an "area 
law", with Qc ^ exp{—cA), with loop area A, and c a 
A'— dependent constant. When K > Kc, Gc decays more 
slowly, only exponentially with the perimeter of the loop. 
What do these two phases correspond to in physical 
terms? Consider first the large K limit, which is the high 
temperature phase of the dual Ising model. As A' ^ 00 
all of the gauge field plaquette sums will be equal to plus 



one. In this case it is possible to choose a gauge in which 
all of the Ising link variables are also unity, Uy = 1. 
In this phase the chargons and spinous can propagate 
at energies above their respective gaps. Apparently, the 
Hamiltonian contains gapped excitations which carry the 
quantum numbers of spinous and chargons. The electron 
has effectively been fractionalized! We denote this exotic 
insulating state with deconfined chargons and spinous as 
X* . It is exceedingly important to emphasize that the 
splintering of the electron into spin and charge carrying 
constituents is conceptually unrelated to the presence or 
absence of spin order. Indeed, electron fractionalization 
can occur even in the presence of strong spin-orbit in- 
teractions which destroys spin-rotational invariance - in 
that case the states of the fermionic /-particles cannot 
be labelled by spin. 

As the coupling K is reduced, so long as the gauge 
theory is in its perimeter phase, the energy to separate 
particles carrying the Z2 charge remains finite, even for 
infinite separation. The chargons and spinous are de- 
confined. Further, with K < cx3, configurations of the 
Z2 gauge theory with plaquette products equal to minus 
one will become possible. One can think of such pla- 
quettes as being "pierced" by non-zero "Z2 - flux" or Z2 
vorticity. Because the number of such plaquettes on any 
given elementary space-time cube is even, the fluxes form 
"tubes" - analogous to Abrikosov vortices in a Type II 
superconductor - which propagate in space-time as par- 
ticles . These particles can scatter and can anihillate in 
pairs, but since their number is conserved modulo 2 they 
carry a conserved Z2 "charge". We will refer to these 
particle-like Z2 vortices as "visons". One can define a 
vison "3-current" , ji, - a field which lives on the links of 
the dual lattice and takes one of two values, zero or one 
- which satisfies, 



(-i)^"=n< 



(91) 



with the plaquette pierced by the dual link. In the de- 
confined phase, I*, these vison particles exist as gapped 
excitations, in addition to the spinous and chargons. In 
terms of the dual Ising model, Sduai, the Ising spins Vi 
are essentially vison creation operators. With the Ising 
model being disordered for large K, the visons (Ising 
spins) are gapped. Thus, the distinct gapped excitations 
in I* are (i) the chargons (ii) the spinous and (iii) the 
visons. An important property of these excitations is 
the existence of long-ranged "statistical" interactions be- 
tween them. Specifically, when a chargon (or a spinon) 
is adiabatically transported around a vison, it acquires 
a geometrical phase factor of tt (because the chargon is 
minimally coupled to the Z2 gauge field). Similarly, a 
vison picks up a tt phase factor upon encircling either a 
chargon or a spinon. Evidently, visons and chargons (or 
spinous) are "relative semions" . 
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As K is reduced further the gauge theory undergoes a 
phase transition at K^ into its "area-law" phase. This im- 
phes that the energy to separate two spinous or chargons, 
inserted as "test" charges at spatial separation, R, grows 
linearly with R. In this "confined" insulating phase, de- 
noted X, free chargons and spinons do not exist in the 
spectrum. The only allowed particle excitations are those 
that are "charge neutral" - that is, invariant under the 
Z2 gauge transformation. Any bound state with an even 
number of chargons plus spinons is "neutral" . In addi- 
tion to the electron, this includes any composite built 
from electrons, such as a Cooper pair or a magnon. In 
the phase I these electron-like excitations will be gapped. 
This phase is the familiar "band insulator" with an even 
number of electrons per unit cell. 

Note that with K < Kc, the dual Ising model orders, 
(vi) ^ 0. This corresponds to a "condensation" of the vi- 
sons. Remarkably, Z2 vortex condensation leads directly 
to a "confinement" for the chargons and spinons. To un- 
derstand confinement directly in terms of the dual Ising 
model, consider the effect of inserting two static "test" 
chargons, separated by a distance R. Each chargon lives 
on a (spatial) plaquette of the dual Ising model. Due 
to the geometrical phase factor between visons and char- 
gons, the presence of a chargon corresponds to a "frus- 
trated plaquette" in the dual Ising model - that is, a pla- 
quette with an odd number of negative Ising couplings. 
To frustrate two plaquettes, it suffices to introduce an 
interconnecting string of negative Ising bonds. In the or- 
dered phase of the dual Ising model, the energy of this 
string will clearly be linear in its length, thereby confin- 
ing the two chargons. 

It is worth drawing a very important distiction between 
the Ising gauge theory considered here, and the gauge 
theories introduced by Baskaran and Anderson ||7| and 
generalized and extensively studied by several authors 
1^. In the simplest version of these theories, the spin 
itself is effectively fractionalizcd, decomposed into a bi- 
linear of spinful (complex) fermion operators, rather than 
splitting the Cooper pair into two chargons as discussed 
above. These spinful fermion operators - the spinons - are 
minimally coupled to a compact U{1) gauge field. But 
in contrast to the Z2 gauge field which exhibits both a 
confined and deconfined phase, the U{1) theory has only 
a single phase |l3]. In this phase, point like monopole 
excitations in 2 -f 1-dimensional space-time always pro- 
liferate, and drive spinon confinement |2q ]. The electron 
is, then, ultimately not expected to be fractionalized in 
these theories. 



2. Superconducting phases 

We now turn to a description of superconductivity 
within the Z2 gauge theory. Since the spinons will be 
gapped into singlets within the superconducting phase, it 



is legitimate to integrate them out, generating once again 
a field strength term for the gauge field as in Eqn. |8^. 
When the dimensonless chargon "hopping" amplitude, 
tc, increases and becomes much larger than unity, one 
expects the chargons to condense, (e"^) ^ 0. For large 
K so that the gauge field is effectively frozen, this char- 
gon condensation transition is simply a 3D classical XY 
transition. Since the chargon carries electric charge e, in 
this phase the charge U{1) summetry is broken, and a 
Meissner effect results. But the chargon also carries Z2 — 
charge, so that the Z2 gauge symmetry is also sponta- 
neously broken. Within a conventional BCS description 
of superconductivity, the order parameter (the Cooper 
pair) carries charge 2e, so one might be tempted to con- 
clude that this "chargon condensate" is perhaps some 
sort of exotic unconventional superconducting phase. In 
particular, it is not a priori clear that the chargon con- 
densate can support a conventional hc/2e BCS vortex. 



To highlight the confusion, it's instructive to focus on 
the regime with large K, where a good description of the 
ground state can be obtained by setting aij = 1 on ev- 
ery link, and taking the chargon phase 4>i a space-time 
independent constant. Consider placing an hc/2e vortex 
at the (spatial) origin. Upon encircling this C/(l) vortex 
at a large distance, the phase of the chargon wavefunc- 
tion must wind by tt. This is of course not possible with 
a smoothly varying phase field, but requires the intro- 
duction of a "cut" running from the vortex to spatial 
infinity across which the phase jumps by tt. The energy 
of this cut is, however, linear in its length with a line 
tension proportional to tc\{e^'^)\'^. It thus appears that 
hc/2e vortices are themselves confined, and not allowed 
in the superconducting chargon condensate. But imagine 
changing the sign of all the Z2 gauge fields, aij, which 
"cross" the cut. This corresponds to placing a Z2 vortex 
at the origin. These sign changes "unfrustrate" the XY 
couplings across the cut, so that the line tension vanishes. 
It is thus apparent that a bound state of a Z2 vortex and 
the hc/2e U{1) vortex (in the phase of the chargon) can 
exist within the chargon condensate. It is this bound 
state which corresponds to the elementary BCS vortex 
in the conventional description of a superconductor. 



It is worth emphasizing that both the "naked" hc/2e 
U{1) vortex and the Z2 vortex - the vison - are confined 
in the superconducting phase. For example, the energy 
cost to pull apart two Z2 vortices also grows linearly with 
separation. To see this, introduce two visons by changing 
the sign of the Z2 gauge field along an interconnecting 
"line" . Due to the chargon condensate which breaks the 
Z2 gauge symmetry making the gauge field "massive" , 
each negative bond costs an energy Ate, implying linear 
confinement. 
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FIG. 2. Schematic zero temperature phase diagram in the 
K — tc plane for local s-wave pairing with an even number of 
electrons per unit cell. 



Thus the distinct massive excitations (apart from the 
Anderson-Higgs plasma mode necessitated by the U{1) 
symmetry breaking) in the chargon condensate are the 
spinons and the BCS hc/2e vortices. This is exactly as 
required in a conventional superconducting phase. Fur- 
ther, since the spinons are minimally coupled to the Z2 
gauge field, there is a long range statistical interaction 
between the spinons and the BCS vortices. In elTect, 
a spinon "sees" the Z2 vortex - which is bound to the 
hc/2e vortex - as a source of "Ising flux". This too is 
as required in a conventional superconductor. Thus, the 
chargon condensate does in fact describe a conventional 
superconducting phase - denoted hereafter as SC. 

A schematic phase diagram is shown in the K—tc plane 
in Fig. ^. The transition from the fractionalized insulator 
X* into SC is essentially a superconductor-insulator tran- 
sition for the charge e chargons. These exist as finite en- 
ergy excitations in I* - superconducting order is obtained 
if they condense. On the other hand, the transition from 
the conventional insulator I into SC can be viewed as a 
superconductor-insulator transition for charge 2e Cooper 
pairs. This can be seen by considering the K = Q limit, 
where it is possible to integrate out the Z2 gauge field 
and arrive at an effective theory of Cooper pair hopping: 



S, 



pair 



-2t2^cos[2{(t)^ 



<!>,)]. 



(92) 



IV. 



ODD NUMBER OF ELECTRONS PER UNIT 
CELL WITH D-WAVE PAIRING 



Having explored the physics of electron fractionaliza- 
tion which follows from the Z2 gauge theory in the sim- 
plest of cases with an even number of particles per site in 



the presence of s-wave pairing correlations, we turn now 
to a much more interesting and challenging situation: 
Correlated Mott insulators with one electron per unit 
cell in the presence of local d-wave pairing correlations. 
As we shall see, in this case the Z2 gauge theory has 
two simple limiting regimes - one describing a d-wave su- 
perconductor and the other a conventional antiferromag- 
netic insulator. But in the interesting crossover regime 
between these two limits, a number of other phases can 
be readily described within the Z2 gauge theory. Besides 
a spin-Peierls ordered phase, the theory gives a simple 
description of the nodal liquid - an exotic fractionalized 
insulator with gapless fermionic quasiparticles. With one 
electron per unit cell, confinement transitions out of the 
d-wave superconductor or nodal liquid are inextricably 
linked to breaking of translational symmetry. 
The full theory of interest can be written as 



S — Sc + Ss + Sb , 

Sc = -2tc 2J 0:ijCOs{(l)i - (j>j), 



(93) 
(94) 



(y) 



As shown in Eqns. O and O, with odd integer A^o there 
is an extra Berry's phase term in the action. 



Sb 



1,3=1- 



(1 



(96) 



It is instructive to consider various limiting cases de- 
scribed by the above action. First consider the limit 
tc = 0. Then Sc = 0, and the (j) fields may be trivially in- 
tegrated out. Surprisingly, the partition function for the 
remaining spin sector of the theory is formally equiva- 
lent to the Heisenberg antiferromagnetic spin model. To 
demonstrate this we first trace over the two allowed val- 
ues of the Z2 gauge field Uij on each link. Consider first 
the spatial links, which enter the action in the form. 



s: 



E 

<.rr'~> T 



/ ^ drr'-^rr' 1 



(97) 



Air' = -KAlrlr' + C.C-) ~ C'ifr^fr'i - (T^i) + C.c). (98) 

For notational simplicity we have suppressed the r index 
on the fermion fields. Tracing over the cr^r' fields for each 
(independent) spatial link and exponentiating the result 
generates a term in the action of the form. 



Ot- — 



^rr'> T 



2_^ In cos/i [Air') 



(99) 



Since A is bi-linear in the fermion fields, upon expand- 
ing in powers of A one generates a series of terms that 
involve multiples of four spinons. 

Now consider the trace of (Jij along the temporal links. 
Recall that the effect of the gauge field aij^i-f along 
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the temporal links is precisely to impose the constraint 
Eqn. pq on the Hilbert space in a Hamiltonian formu- 
lation. With the (j) fields integrated out, at tc = 0, this 
constraint reduces to requiring 



(-1)"/=-!, 



(100) 



at each site of the spatial lattice. Due to Pauli exclusion 
this is equivalent to the constraint that n/ = 1 at each 
site. Thus, after tracing out the a field, the Hamiltonian 
obtained from Sr is constrained to operate on a Hilbert 
space with exactly one spinon per site. This Hamiltonian 
consists of a sum of terms for each nearest neighbour spa- 
tial link. With the additional requirement of spin rota- 
tion symmetry, the Hamiltonian must take the form of 
the Heisenberg spin Hamiltonian, 



11 ij y ij f ' ij r^' 



(101) 



<rr'> 



This can be verified directly from Sr by expanding out 
the Incosh term, and re-expressing the spinon operators 
in terms of the spin operators, 5'^ = f^crfr- This leads 
to an explicit expression for the exchange interaction: 



..!(,<.,. 2^ 



(102) 



where e is the discrete time slice defined in Eqn. W4. 

Recovering the Heisenberg antiferromagnet in the limit 
tc — > provides a way to obtain a rough estimate for the 
saddle-point parameter xo- First, we note that i* and t 
can be re-expressed in terms of the parameters t, u, J, A 
and Xo using Eqns. [79|and|80|. Though these relations are 
strictly valid for s-wave pairing, they suffice to give rough 
estimates even for the d-wave case. It is, however, nec- 
essary to modify the equation for t^ due to the slightly 
different decoupling in the d-wave case (See Appendix 
y|). Assuming that the saddle point value tjq ~ xo; "^e 
get 



t' 



A 



rt". 



(103) 



Combining Eqns. |79| and |103| with Eqn. |102| and assum- 
ing A << J, leads to an estimate for xo^ 



Xo^ \-pT 



(104) 



which is appropriate in the limit of large u/t. Having es- 
timated xoi one can use Eqns. fF 



and 
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to obtain 



estimates for the three dimensionless coupling constants, 
tj,, V^ and t^, respectively. The resulting estimates are 
given in Eqn. EG. 

Having established the equivalence of the action in 
Eqn. p3 to the Heisenberg antiferromagnet in the limit 
tc —^ 0, we briefly consider the opposite large tc limit. 



With sufRciently large tc the chargons will condense, and 
as argued in the previous section this describes a conven- 
tional superconducting phase. But due to the assumed 
form of the pairing correlations, the pairing symmetry 
here will be d^2_y2 . Thus, the Z2 gauge theory in Eqn. B3 
has the remarkable property that it describes a conven- 
tional antiferromagnet for small chargon coupling, and a 
conventional d^^^y^ superconductor in the opposite ex- 
treme. We now turn our attention to the exceedingly 
interesting regime between these two limits. 



3. Correlated Mott insulators 

When the chargon coupling strength tc is small, the 
chargons will be gapped out, and the system in an insu- 
lating phase. In this case, it is appropriate to integrate 
out the chargon fields to obtain an effective action for the 
spinous and the gauge field a. The main result of this 
integration will be to generate a plaquette product term 
of the form, 



S. 



-^Etn-. 



(105) 



The full remaining action which is valid within the insu- 
lating phases is then simply. 



S — Ss -\- Sa + Si 



(106) 



The parameter K depends on the coupling tc, vanish- 
ing at tc = and increasing monotonically with tc- The 
transition to superconductivity will occur when tc ^ 1. 
Near this limit, but on the insulating side, the value of 
K will also be of order one. Keeping this in mind, we 
first find it convenient to analyze the phase diagram of 
the above action for arbitrary K, incorporating later the 
superconducting phase. 
The action in Eqn. 
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has three dimensionless cou- 
pling constants, t^,t^ and K. Considerable progress can 
be made in determining the phase diagram by focussing 
on three different limits. The first, considered above, is 
K ^ where the model reduces to the Heisenberg spin 
model. The second tractable limit is large K. li K = 00 
the gauge field is frozen out and it is possible to choose 
a gauge with aij = 1 on every link. Then, the only 
remaining piece of the action describes non-interacting 
spinous with a gapless "d-wave" dispersion at four points 
in the Brillouin zone. This is the "nodal liquid" phase - 
obtained in earlier work ^^M by vortex-pairing within 
a dual vortex formulation. The nodal liquid is a frac- 
tionalized insulator with deconfined, gapless spinous and 
gapped chargons. For large but finite K and in the ab- 
sence of Sb, the Z2 gauge theory is in its perimeter law 
phase. As we show below, this continues to hold even in 
the presence of Sb - in fact, the region of stability of the 
perimeter phase is enhanced by the Sb term. Thus, the 
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chargons and spinons remain deconfined and the nodal 
liquid phase survives for large but finite K. 

As with the fractionalized insulator discussed in See- 
apart from the chargons and the spinons there 
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tion 

is another distinct excitation in the nodal liquid phase 
- the Z2 vortex configuration in the a field, dubbed the 
"vison" . The vison is a gapped excitation in the nodal 
liquid. As before, due to the minimal coupling of the 
chargons and the spinons to the Z2 gauge field a, they 
each acquire a phase of tt upon encircling a vison. There 
is thus a long ranged statistical interaction between a 
chargon (or a spinon) and a vison. 

The third tractable limit of the action Eqn. 
small t^ and t^ . (Estimates appropriate to the cuprates 
obtained from Eqn. g^, suggest that these couplings will 
most likely be much smaller than one.) In the extreme 
limit oit^ = t^ = 0, we are left with a pure Z2 gauge the- 
ory described by Seff — Sa + Sb- To explore the effects 
of the Berry's phase term Sb on the gauge theory, it is 
useful to pass to the dual representation. Recall that for 
5_B = the dual theory is simply the 2 -I- 1— dimensional 
Ising model, with the Ising spin operators {vi = ±1) cre- 
ating the "vison" excitations. To implement the duality 
transformation with the Berry's phase term present, it is 
convenient to first rewrite it as. 






;)(i-n^^f)- 



(107) 



Here /z?^* can be viewed as an "external" Z2 gauge field 
living on the links of the dual lattice, which satisfies 
ria A*?/ * = — 1 through every spatial plaquette. In this 
form one can readily generalize the duality transforma- 
tion in Appendix O to give. 



Sdu 
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with dual coupling satisfying; tanh{Kd) — e^^^ . Due 
to the Berry's phase term, every spatial plaquette (with 
normals along the time direction) in the dual Ising model 
is frustrated. In the time continuum limit this becomes 
a 2d quantum transverse-field Ising model which is fully 
frustrated. 

The quantum Ising model on a fully frustrated square 
lattice has been studied extensively by several authors 
|g9|,|6[. In particular, Jalabert and Sachdev |1^ studied 
the model numerically (not coincidentally) in the con- 
text of frustrated quantum Heisenberg spin models. For 
small Kd the Ising model exhibits the usual paramagnetic 
phase, in which the visons are gapped (uncondensed) 
with (uj) = 0. This corresponds to the "low tempera- 
ture" phase of the gauge theory. Deep within this phase 
one can set aij = 1 on all the links, which implies (for 
t'',i 7^ 0) that the chargons and spinons are deconfined. 
This is the nodal liquid phase discussed earlier. It is note- 
worthy that the frustration in the Ising model - which is 



a direct consequence of being in a Mott insulator with 
one electron per site - enhances the stability of the frac- 
tionalized nodal liquid phase (the paramagnetic phase of 
the Ising model). 



FIG. 3. One possible ordered low temperature phase of the 
fully frustrated transverse field Ising model in two spatial di- 
mensions. The thick lines represent the frustrated bonds. The 
dashed lines denote the links of the dual lattice where the cor- 
responding "singlet bonds" live. 



As Kd is increased, it has been found ||T3] that the Ising 
model orders - breaking the global Z2 spin flip symmetry. 
But due to the frustration, this ordering is accompanied 
by a spontaneous breaking of translational symmetry. It 
is convenient to characterize this symmetry breaking in 
terms of the gauge- invariant energy densities of the near- 
neighbor bonds: £ij = —{vifj,f^*Vj). It is found that some 
of the bonds are "frustrated" with positive £y , while the 
remaining are "happy" with negative bond energies. In 
the spatially broken ordered phases, it is fomid that these 
frustrated bonds form lines (see Fig. ^), which run along 
the principle axis of the square lattice (columns or rows). 
There are four favored configurations, corresponding to 
frustrated bonds along every other column, or along ev- 
ery other row. Within each of these phases, a particular 
gauge choice can be made with /i^^' = — 1 on each "frus- 
trated" bond. With this choice of gauge, the Ising spins, 
Vi, exhibit a global ferromagnetic ordering. Altogether, 
the ground state is thus eight-fold degenerate and breaks 
the Z2 spin flip, translational and rotational symmetries. 

In general, several other ordered phases of the fully 
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frustrated Ising model are possible - some of these are 
explored in the Landau theory of the first reference in 
Ref. Ig^. These phases may perhaps be stabilized by 
very large Kd, and/or longer ranged interactions beyond 
the simplest nearest neighbour model studied in Ref. [^ . 
We will not consider these other possibile phases in the 
present paper. 

What are the effects of a small non-zero i" and t^ 
which couple the spinous to the Z2 gauge field? In the 
context of quantum antifcrromagnets, Sachdev [nqJlTl 
has suggested that the spatial ordering of the Ising model 
corresponds to a spin-Peierls ordering. This interpreta- 
tion appears to be consistent within our present frame- 
work. Specifically, associated with each frustrated bond 
in the Ising model, is a corresponding frustrated plaque- 
tte on the dual lattice "pierced" by that bond. The ex- 
pectation value of the plaquette product in the gauge the- 
ory will therefore be modulated in these ordered phases. 



with (Hn <yv 



~£ij- Upon including the coupling to 



the spinous, this modulation of the energy density will, 
in general, induce a modulation in various other physical 
quantities. In particular, the quantum expectation value 
< Sr ■ S,' > evaluated for each bond will be spatially 
modulated - bonds which "cross" the frustrated lines of 
the dual lattice will have a different value for this expec- 
tation value from other bonds. Presuming the spin rota- 
tion invariance remains unbroken, this state corresponds 
to a spin-Peierls phase - which we denote as SV. The 
"singlet bonds" in this phase are arranged in a columnar 
fashion - running perpendicular to the lines of frustrated 
bonds in the dual Ising model as depicted in Fig. 0. 

Since the Ising spins in the fully frustrated Ising model 
order ferromagnetically in these modulated phases (with 
an appropriate gauge choice for /if,^') implying a vison 
condensation, (vi) 7^ 0, confinement is expected. To see 
this, consider evaluating the Wilson loop correlator de- 
fined in Eqn. ^. In the dual frustrated Ising model, this 
corresponds to changing the sign of all the Ising couplings 
on bonds which "pierce" throught the loop. Being fer- 
romagnetically ordered, this will cost an energy (action) 
proportional to the area of the loop - the signature of 
confinement. Thus, as expected, the spin-Peierls state is 
a conventional insulator, with confined spinous and char- 
gons. The gapped spin-one excitations made by breaking 
the singlet bonds can then be thought of as a (confined) 
pair of spinous. 

The three limiting cases discussed above suggest the 
phase diagram shown in Fig. H for the action in Eqn. 106 
Consider first the regime with small t'' and t^. At very 
small K a conventional antiferromagnetic insulator is ex- 
pected. With increasing K there is presumably a phase 
transition into a conventional spin-Peierls insulator with 
confined chargons and spinous. Upon further increas- 
ing K, the spin Peierls phase undergoes a deconfinement 
transition into the fractionalized nodal liquid phase. For 
large t'' and f^, the antiferromagnet and nodal liquid 



phases will still be present in the limit of very small and 
large K, respectively. But it is not clear which phases 
will be present when all three of the coupling constants 
are of order one. In particular, it is unclear whether it 
is possible to have a direct second order phase transi- 
tion from the antiferromagnet into the nodal liquid , or 
whether there will always be an intervening (spin-Peierls) 
phase. 
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FIG. 4. Schematic zero temperature phase diagram show- 
ing one possible scenario for the evolution from the anti- 
ferromagnet (AF) to the d-wave superconductor dSC. In 
this scenario, all the insulating phases are conventional. The 
thick lines indicate confinement of the chargons and spinous. 
For concreteness, we have chosen to display a particular se- 
quence of confined phases, namely, a transition from AF to a 
spin-Peierls (SP) insulator, and a further transition to dSC. 
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FIG. 5. The other qualitatively different scenario for the 
evolution from the antiferromagnet to the d-wave supercon- 
ductor. In this case, on increasing tc, a transition to the frac- 
tionalized nodal liquid (NL) phase occurs before the onset of 
superconductivity. 



We now discuss the implications of these results for 
the phase diagram of the full Z2 gauge theory in which 
the charge degrees of freedom are present and supercon- 
ductivity is possible. Of primary interest is the evolution 
from the antiferromagnet to the d— wave superconductor 
upon increasing the chargon coupling, tc- A transition 
into the superconductor is expected to occur at some 
critical chargon coupling, t*, of order one. For smaller 
tc in the insulating regime, the dimensionless coupling 
K will at most be of order one. One can imagine two 
qualitatively distinct possibilities upon tuning towards 
the superconductor from the insulating phases. First, it 
may be that even when tc increases to t* , the value of K 
will remain smaller than the critical value needed for de- 
confinement, Kc- In this case, all the intermediate phases 
between the antiferromagnet and the superconductor will 
be conventional confined phases. This is illustrated in 
Fig. ^. Alternately, it may be that K exceeds Kc before 
the onset of superconductivity. This would imply the ex- 
istence of the deconfined nodal liquid phase intervening 
between the d— wave superconductor and a conventional 
insulator. This is illustrated in Fig. 0. 

Which one of these two possibilities is realized will pre- 
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sumably depend sensitively on microscopic details. In- 
deed, since K is of order one when tc approaches t*, it 
seems likely that the onset of superconductivity will oc- 
cur close to the boundary between the confined and de- 
confined insulating phases. But in any event, our analy- 
sis has firmly established the possibility of the deconfined 
nodal liquid phase. It remains as a challenge to determine 
whether this exotic fractionalized insulator is realized in 
the cuprates. 

In the next section, we will describe much of the 
physics discussed here in a dual formalism in terms of 
vortices rather than the chargons. This will provide con- 
siderable further insight, and make connections with ear- 
lier approaches. 



V. DUAL VORTEX REPRESENTATION 

For a system of interacting bosons in two spatial di- 
mensions, it is well-known that the insulating phases can 
be described as a condensate of vortices in the many par- 
ticle wave-function. More formally, it is possible to set 
up a dual description where the vortices, rather than the 
particles, are the fundamental degrees of freedom such 
that the insulating phase is a vortex condensate while the 
superfluid phase is the vortex vacuum. For the electronic 
systems considered in this paper, it is natural to attempt 
to do the same, and work with a dual description in terms 
of vortices in the Cooper pair phase (p, and the spinous. 
Since the Cooper pair has charge 2e, these are the ^ 
vortices which occur in a conventional superconductor. 
Besides providing additional insight into the mechanism 
and nature of electron fractionalization, passing to a dual 
vortex description enables us to make contact with earlier 
work which describes fractionalization in terms of vortex 
pairing. 

We will start with the full chargon-spinon action S — 
Sc + Ss + Sb discussed in the last section, and perform a 
duality transformation to trade the chargon fields for the 
■|^ vortices. This differs somewhat from the conventional 
duality transformation [H from bosons to vortices due 
to the coupling of the chargons to the Z2 gauge field. 

To understand how to deal with the chargon coupling 
to the a field, it is useful to first review the well-known 
self-duality of the Z2 gauge theory with Ising matter 
fields in 2 + 1 dimensions. This is done in detail in Ap- 
pendix p. The duality proceeds by first rewriting the 
partition function in terms of a Z2 current for the Ising 
matter fields and the Z2 gauge field, Uij. The Z2 cur- 
rent lives on the links of the lattice and can take one of 
two values 0, 1. It is conserved modulo 2 at each site of 
the lattice. This conservation law can be implemented 
by writing the Z2 current as the flux of a dual Z2 gauge 
field, denoted as /Xy. (This is completely analogous to 
the duality of the three dimensional classical XY model) . 
Eliminating the Z2 current in favor of the dual gauge field 



gives an action written entirely in terms of two Z2 gauge 
fields {aij and /iy) which are duals of each other. The 
original Z2 gauge field, cfij may be eliminated by express- 
ing its flux as the current of a dual Ising matter field, the 
vison Vi. The resulting theory has the same form as the 
original Z2 gauge theory with matter fields, but is dual 
to it. 

To obtain a dual representation of the system of char- 
gons and spinous coupled to the Z2 gauge field aij, we 
need to combine the dual representation of the Z2 gauge 
theory with the standard duality transformation of the 
XY model. As shown in detail in Appendix H this is 
readily done. For the time being, we will only consider 
the situation with local d-wave pairing and an odd num- 
ber of electrons per unit cell. The result is a lattice action 
in terms of ^ vortices, which are minimally coupled to 
a fluctuating C/(l) gauge fleld a whose circulation is the 
total electrical current. In addition, the hc/2e vortices 
are minimally coupled to a Z2 gauge field ^ij. The full 
action is given by 



S — Si, + Sa + Ss + Scs + Sb , 
Sv = -ty ^ ^ijCos{6i - 6j + -^), 



(109) 
(110) 
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^'^ = ^E(^x «■'.)'' (Ill) 
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^<^^=E4(i-n^)(i-/^^^) 



(113) 



Here e**' creates the || vortex, and fi is the spinon as 
before. The first term represents single vortex hopping, 
while the second is a kinetic term for the U{1) gauge 
field aij . The fiux of a is the total electrical current - in 
particular a flux of 2tt through a spatial plaquette adds 
an electric charge of one - a chargon. Together these 
two terms comprise the usual dual vortex representation 
of a set of charge 2e Cooper pairs, except that here the 
vortices are minimally coupled to an additional Z2 gauge 
field /iij. This leads to a vortex-spinon coupling medi- 
ated by Scs- This term has a structure very similar to 
a Chern-Simons term (although it is for the group Z2), 
and as discussed below, plays a similar role. The Berry's 
phase term Sb is the same as before. 

The full dual action is invariant under a local U{1) 
gauge transformation 



0, + A,: 



A,; - A, 



(114) 
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This is standard in the dual vortex description of XY 
models in three dimensions. The corresponding con- 
served charge is the vorticity The action has an ad- 
ditional Z2 gauge symmetry under which 
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jllj > Ciflij€j. 



(116) 



A. Phases 



with ej = ±1. We emphasize that this gauge symme- 
try is distinct from the local Z2 gauge symmetry of the 
spinon-chargon action, but in fact is dual to it. 

To get some intuition for the term Scs^ it is instruc- 
tive to replace the vortex hopping term in the action by 
a Villain potential, 



^t^cosOij 



-Jl/2t^ ij^etj 



J„ = -c 



(117) 



where Qij = 6i — 



'j + ^ -t- f (1 - ^y) is the gauge in- 
variant phase difference. Here the integer field Jy that 
lives on the links of the lattice represents the 3-current of 
the hc/2e vortices. After this replacement it is possible 
to explicitly perform the summation over the gauge field 
fj,ij. For each link of the lattice this contributes a term 
to the partition function of the form, 1 -I- (— l)"^" Yla cr, 
which vanishes unless 



(-i/-n- 



(118) 



Thus, the Chern-Simons term has effectively attached a 
Z2 flux of the gauge field a - a "vison" - to each hc/2e 
vortex. As discussed in Section II], this composite com- 



prised of an hc/2e vortex bound to the Z2 vison is nothing 
but the familiar BCS vortex. Due to the attached vison, 
when a spinon is taken around the BCS vortex it acquires 
the expected n phase factor. 

Alternatively, it is possible to perform an "integration 
by parts" on Scs which effectively exchanges the role of 
a and /i, and then perform a summation over a. This 
leads to the additional constraint. 
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with Jf the spinon 3-current. A Z2 flux in the gauge field 
fj, has thereby been attached to each spinon. More pre- 
cisely, since the spinon number is only conserved modulo 
2 due to the anomalous pairing term, the Z2 flux is at- 
tached whenever an odd number of spinous propagates. 
The net effect of this Z2 Chern-Simons term is to imple- 
ment mathematically the long-ranged statistical interac- 
tion between BCS vortices and spinous. This kind of flux 
attachment may be familiar to many readers for the U{1) 
group from theories of the quantum Hall effect. But since 
the spinon number itself is not conserved, implement- 
ing this statistical interaction with a C/(l) Chern-Simons 
term is problematic. It is a remarkable aspect of the du- 
ality transformation in Appendix O, that this Ising-like 
Chern-Simons terms emerges so naturally. 



We now analyze the phases in this dual vortex descrip- 
tion, focussing on the most interesting case of an odd 
number of electrons per site with local d-wave pairing 
correlations. In the vortex description the superconduct- 
ing phase corresponds to a vortex vacuum, and the insu- 
lating phases are vortex condensates. We consider flrst 
two simple limiting cases, firstly the superconductor with 
vanishingly small vortex hopping <„ —f 0, and then the 
insulator with t„ — > 00. 

When ty is zero the summation over the gauge field 
fi can be performed, giving the constraint Yin c = 1. 
It is then possible to pick a gauge with cr^ = 1 on ev- 
ery link. The resulting action has two pieces, Sa which 
describes the gapless sound mode of the superconduc- 
tor (gapped when long-ranged Coulomb interactions are 
included) and the spinon piece Sg- With aij = 1 the 
spinous can freely propagate and describe the gapless 
nodal quasiparticles. A correct description of a conven- 
tional d-wave superconductor is thereby recovered. 

Consider next the opposite limit with i„ — > cxd. In this 
regime the hc/2e vortices will condense, (e**') 7^ 0. The 
dual "Anderson-Higgs" mechanism leads to a mass term 
for the gauge field a^ , indicative of a charge gap. With 
one electron per unit cell the resulting phase is thus a 
Mott insulator. In the absence of any gapped charge ex- 
citations (Ax a = 0), it is possible to choose a gauge with 
ttij = on every link. The vortex hopping term becomes, 
Sv = —hJ2ij t^ij^ with a non-zero "field": h = it,|(e'®')p. 
When this "field" is large one can set /iy = 1 on each 
link, so that the Chern-Simons terms vanishes. The full 
action then reduces to Seff — Ss + Sb- At this stage the 
summation over the a gauge field can be performed ex- 
plicitly. As detailed in the previous section, the resulting 
model reduces to a simple 2d near-neighbor Heisenberg 
antiferromagnet. Thus, we readily recover the simple an- 
tiferromagnet from the dual representation by condens- 
ing hc/2e vortices. 

Finally, in this section we wish to recover a dual de- 
scription of the fractionalized "nodal liquid" . Since the 
nodal liquid is electrically insulating it requires vortex 
condensation. But as established in the previous section, 
the nodal liquid supports gapped Z2 vortices - the "vison" 
excitations. Since the Chern Simons term attaches a vi- 
son to each hc/2e vortex, it is clear that to obtain the 
nodal liquid the hc/2e vortices cannot be condensed. But 
since the square of the vison operator is unity {vf = 1), a 
pair of hc/2e BCS vortices does not carry a vison with it. 
As we now show, the nodal liquid can be obtained from 
the d-wave superconductor by pairing BCS vortices, and 
then condensing the hc/e vortex composite. 

To this end, we add an extra vortex pair hopping term 
to the action, 
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Here, e*^^' = (e'^')^, thus creating a pair of BCS vortices. 
Notice that the hc/e vortex is also minimally coupled to 
the [/(I) gauge field - as required by the dual U{1) sym- 
metry of the action - but is not coupled to the Z2 gauge 
field, fj,ij, because it carries no vison charge. We now 
consider taking t2v large and condensing the hc/e vor- 
tex, (e*^^') ^ 0, keeping the hc/2e vortex uncondensed. 
Before doing this it is convenient to re-express the hc/2e 
vortex as 



pie2,/2 



(121) 



with Vi — ±1 the vison operator. Notice that with this 
identification the field 62 can be treated as an angular 
variable, since the right side is invariant under the com- 
bined transformation, 82—^02 + 27r and Vi — )■ — Wi. We 
finally find it convenient to absorb the field d2i into the 
gauge field aij by the gauge transformation. 
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In this gauge, the vortex hopping terms become 

Sv = -ty ^ VifiijVjCos (-^j , (123) 

S2v = -hv ^ cos{aij). (124) 

In the insulating phase with large t2v there will again 
be a charge gap due to the dual Anderson-Higgs mech- 
anism, coming from the hc/e vortex condensate. Above 
the gap will be charge e chargons, corresponding to a 27r 
fiux tube in a^. In the absence of any charged excita- 
tions one can set aij = 0, and the single vortex hopping 
term becomes, 



Oy 



Zy y ^ VifiijVj . 
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The full effective action is Seff = Su + Sg + Scs + Sb- 
When ty is small the visons will be uncondensed (vi) — 0. 
In this limit the summation over the fi gauge field can 
be performed, and due to the Chern Simons term leads 
to the constraint, Jin ''' — 1- One can then choose a 
gauge with aij = 1 on each link, which sets Sb — 0. The 
only remaining term in 5*6// describes free propagating 
spinous. These are the gapless nodons in the insulating 
nodal liquid. 

We thereby recover a decription of the nodal liquid 
from the dual vortex formulation. In addition to the gap- 
less nodons, the nodal liquid supports two gapped excita- 
tions; the chargon and the vison. As clear from the above 
analysis, the vison is simply a remnant of the hc/2e BCS 
vortex which survives into the nodal liquid upon con- 
densation of the hc/e vortex pair. Physically, since the 



vorticity is only conserved modulo 2 (in units of hc/2e) 
once the field e*^^' has condensed, only a conserved Z2 
remains from the hc/2e BCS vortex. As before, the vi- 
son picks up a TT phase change when it is transported 
around either a spinon or a chargon. To see this, note 
that a chargon corresponds to a tt flux in aij/ 2 and the 
nodon (spinon) a n fiux in /x^ . As seen in Eqn. 



122, the 



vison is minimally coupled to both of these gauge fields, 
thus acquiring a sign change upon encircling the spinon 
or chargon. 

It is worth emphasizing that a clear mechanism for vor- 
tex pairing can be found from the analysis in the previous 
section. Since the chargons and visons (or vortices) have 
a long-ranged statistical interaction, motion of the charge 
is greatly impeded by the presence of unpaired visons. 
On the other hand, once the hc/2e vortices are paired, 
the charge can move coherently. Thus, the presence of 
a large kinetic energy makes vortex pairing energetically 
favorable. 

It is finally worth mentioning that in the limit Ss — 0, 
one readily recovers the fully frustrated Ising model con- 
sidered in the previous section. To see this, note first that 
Sb can be re- written in the form of a Chern-Simons terms 
with n replaced by /i'^^*, where Hn A^*^"^* — ^1 through 
all spatial plaquettes. With Sg — 0, one can then per- 
form the summation over the a gauge field, and this sets 
f-ij = MiT'- The remaining term in Seff is the fully frus- 
rtated Ising model: 
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VI. DOPING 

Our analysis has so far focussed only on situations with 
an integer number, N^^ of electrons per unit cell. Finite 
doping leading to non-integer Nq does not crucially mod- 
ify our discussion of fractionalization issues. Indeed, both 
confined and fractionalized insulating phases can exist for 
non-zero doping. At a qualitative level, in both kinds of 
insulating phases, the main effect of non-integer A'o will 
be to induce charge order, accompanied by translational 
symmetry breaking. The precise nature of this charge or- 
der presumably depends on the details of the system, and 
may be sensitive to the presence of long-ranged Coulomb 
interactions. 

Formally, non-integer values of Nq can be incorporated 
into either the particle or vortex representations as fol- 
lows. In the particle representation, as discussed in Sec- 
tion n, the main effect of non-integer TVq is to modify the 
Berry phase term to 



Sb — — 






No{2Trhj--il-a,,)). 



(127) 



20 



Here, hj is an integer defined on each temporal link given 

by 



h 



Int 



27r 
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where $ij = (j)i — (j)j + ^(1 — aij) is the gauge-invariant 
phase difference between two sites. When A^o is not an 
integer, this Berry phase term leads to complex Boltz- 
mann weights in the partition function sum. This is not 
too surprising- even in the absence of any gauge field cou- 
pling, the partition function for simple Boson Hubbard 
models at arbitrary chemical potential involves complex 
weights. 

The presence of such complex weights does not pose 
a problem for the existence of the fractionalized insula- 
tor. We recall that the fractionalized phase is obtained 
when the gauge field aij is in its perimeter phase. Deep 
in this phase, we may set aij « 1 on each space-time 
link so that the Berry phase term Sb becomes indepen- 
dent of aij . The resulting action then describes a lattice 
model of bosonic chargons at filling A^o a-nd the fermionic 
spinous, decoupled from one another. Thus the char- 
gons and spinous will still be deconfined. However, the 
ground state will generally exhibit charge ordering ac- 
companied by broken translational invariance. Confined 
conventional insulating phases at non-integer A'o clearly 
also exist. 

Numerical simulations of the Z2 gauge theory at arbi- 
trary Nq to determine the precise nature of the charge 
ordering in these insulating phases will be seriously ham- 
pered by the presence of these complex weights in the 
partition function. Fortunately, in the dual vortex rep- 
resentation, non-integer Nq enters in a more innocuous 
manner. To generalize the duality transformation to ar- 
bitrary A^o is straightforward, because the Villain rep- 
resentation of the chargon hopping term in Eqn. 
simply modified to read, 
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Here, iV^ = A'o for temporal links, and is zero other- 
wise. Proceeding with the duality transformation gives 
the action. 



For the particular case of odd integer A^O: it is instruc- 
tive to see how the term Sb may be recovered. To that 
end, we define a new "external" gauge field o"^^* on the 
links of the dual lattice such that, 



Ax a? 



27riVi, 



(132) 



We now absorb a*^^* into a by the shift a ^ a ~ a*^^*. 
This eliminates a*^^* so that Sa — > Sa, but modifies the 
vortex hopping term which becomes. 
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For odd integer A'o (say A'o = 1) one may choose, 



</* =27rny. 



(134) 



with integer riij , which satisfies A x Uij = Ao = 1 for ev- 
ery spatial plaquette and is zero for all other plaquettes. 
With this choice we may write. 



^v — f^v / f^ij t^ij eOS[ 



'' + T>' 
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where /xfj* = (— 1)"'j. Notice that the fiux Hn M is 
— 1 for every spatial plaquette, and zero for other pla- 
quettes. If we now perform the shift, fj, —^ /i^*^^*, 
the field /i*^^* is eliminated from Sy but reappears in 
S'cs(/i/i^^*). But upon noting the form of the Berry's 



phase term in Eqn. 107, one can easily demonstrate that 
5cs'(/^m'^^*) = 'S'cs'(m) + Sb- We thereby recover the ear- 
lier Berry's phase form for the case with odd integer A'o. 
The dual representation for arbitrary A'o is simpler 
looking than the one in the particle formulation, and is 
probably better suited to discuss issues such as the na- 
ture of charge ordering at finite doping. In particular, if 
we ignore the coupling to the spinons and set Od M = 1' 
the remaining partition function sum involves only real 
weights, and can presumably be evaluated numerically. 



S^Sy + Ss + Scs + Sa: (130) 

where the first three terms are the same as before in 



Eqn. 109 The last term, which was equal to Sa + Sb for 
integer Aq, becomes instead. 
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Notice that in this dual representation, Aq acts like an 
external "magnetic field" piercing each spatial plaquette. 



VII. OTHER EXOTIC FRACTIONALIZED 

PHASES 



In this section we will briefiy explore the possibility of 
obtaining other fractionalized phases different than the 
ones discussed so far. The most interesting phase that 
emerges is a novel fractionalized superconductor - we will 
describe its properties in both the particle and vortex 
formulations. 
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A. Particle description 

In earlier sections we argued that when the charge e 
chargons condense, the resulting phase is a conventional 
superconductor. This is perhaps surprising, since in a 
conventional BCS description the order parameter car- 
ries charge 2e. One might ask whether it is possible to 
have a superconducting phase in which the chargon pairs 
{i.e the Cooper pairs) have condensed, while single char- 
gons have not. As we now demonstrate, such a super- 
conducting phase - which we denote as SC* - can exist 
and has a surprisingly simple description in terms of our 
Z2 gauge theory. For simplicity, we will initially present 
the discussion for s-wave pairing with an even number of 
electrons per unit cell. 

The appropriate action from Eqn's. ^ and [Sg in Sec- 
tion III, takes the form S = Sc + Ss + Sk- As discussed 
there, the kinetic term for the gauge field 5^, although 
not present in the original action, will in any case be 
generated upon integrating out high-energy modes. To 
access the chargon pair condensate phase, it is extremely 
convenient to add an explicit pair hopping term to the 
action, Spair from Eqn. 93. For large pair-hopping am- 
plitude, ^2, the chargon pairs will condense, leaving the 
single chargons uncondensed: 
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This still breaks the global U{1) charge symmetry, and 
so describes a superconductor, but one with rather ex- 
otic properties. To examine this phase it suffices to take 
^2 ^ 00 which allows one to set 20^ equal to 27r times an 
integer, or equivalently. 
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with the value of the Ising spins, Si = ±1, arbitrary. 
In this limit, the chargon creation operator equals the 



Ising spin: e 



i<f)i 



Si. After integrating out the massive 



spinous, this leaves an effective theory of the form: 
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with tc the chargon "hopping" strength. 

This theory, which describes Ising spins "minimally 
coupled" to a Z2 gauge field, has been extensively studied 
by Fradkin and Shenker pfl] a-s a toy model of confine- 
ment. The phase diagram in the tc — K plane is shown in 
Fig. 0. In the K ^ 00 limit the model reduces to a global 
Ising model for the spins. With increasing tc there is an 
Ising transition into a phase with (sj) 7^ (the "Higgs" 
phase) , which corresponds to the chargon-condensed SC 
phase. Along the tc = axis the pure Z2 gauge field ex- 
hibits a confinement transition with decreasing K. Frad- 
kin and Shenker argued that the "Higgs" and confined 
phases could be continuously connected, by noting the 



absence of a phase transition along the tc = 00 and K — 
lines. Moreover, as detailed in Appendix y, this model is 
in fact self-dual, and maps into an equivalent model with 
new parameters reflected across the dashed line. 



00 








FIG. 6. Schematic zero temperature phase diagram for the 
Z2 gauge t heory coupled to matter fields described by the 
action Eqn. 138 . 



The phase with large K but small tc corresponds to the 
exotic new superconducting phase, SC* . In this phase 
there are four deconfined massive excitations: (i) the 
spinon, (ii) an hc/2e U{1) vortex, (iii) the Ising spin Si 
and (iv) the Z2 vortex in the gauge field a - the "vison" . 
In striking constrast to a conventional superconducting 
phase, in SC* the U{1) and Z2 vortices can exist as sep- 
arate excitations, and are not confined to one another. 
In order to distinguish this hc/2e vortex from the BCS 
vortex, we will refer to it as an hc/2e vorton. The Ising 
spin excitation s is a remnant of the chargon. In the 
paired-chargon condensate SC* phase, the global U{1) 
charge symmetry is not fully broken - there is an unbro- 
ken Z2 "charge" symmetry (s, — > — s,) corresponding to 
an invariance under a sign change of the chargon oper- 
ator. Although the electrical U{1) charge of the char- 
gon is not conserved, the chargon number is conserved 
modulo 2, a refiection of this unbroken Ising symme- 
try. Indeed, one can define a conserved Ising "charge" 
as, Q2 — (—1)^ — ±1, where N is the chargon num- 
ber operator. Since the Ising spin operator changes the 
sign of Q2, the massive spin excitation carries the con- 
served Z2 electrical charge of the chargon. We refer to 
this excitation as an "ison" . 

To gain some physical insight into this strange ison par- 
ticle, consider what happens when an electron is added to 
a superconductor. The electron creation operator can be 
decomposed into the product of a spinon and a chargon. 
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The second equality is valid within the two supercon- 
ducting phases. In the conventional superconductor SC, 
the ison is also condensed, (si) ^ 0, so that the electron 
is essentially equal to the spinon. Thus the spin of the 
added electron is carried away by the spinon - the conven- 
tional BCS quasiparticle - whereas the electrical charge 
is carried by the condensate. On the other hand, in SC* 
adding an electron not only increases the conserved spin 
by 1/2, but changes the conserved Z2 "electrical charge" . 
The spin and Z2 charge are carried away by two separate 
massive excitations - the spinon and ison. Thus, the SC* 
phase exhibits an exotic form of spin-charge separation. 

It is again important to ask about geometric phase fac- 
tors acquired when any of the four massive excitations in 
SC* encircle another. First, note that both the ison and 
the spinon are minimally coupled to the gauge field a. 
Consequently, they both acquire a phase factor of n on 
encircling the Z2 vortex, namely the vison. The ison, 
being a remnant of a chargon, also acquires a phase of tt 
on encircling an hc/2e vorton. Thus the pairs - (spinon, 
vison), (ison, vison), and (ison, hc/2e vorton)- acquire 
phase factors of tt upon encircling one another. Equiv- 
alently, there are long-ranged statistical interactions be- 
tween any two members of a pair. All other pairs of 
excitations do not acquire any geometrical phase factors. 
Note in particular that the hc/2e vorton, being unbound 
from the Z2 vison, does not have a long range statistical 
interaction with the spinon in SC* . This distinguishing 
feature will have several important consequences in the 
dual vortex description developed in the next section. 



Z2 vortex). The result is the BCS hc/2e vortex, as dis- 
cussed earlier in Section III. 



SC 



SC 



The transition from I* into SC* upon increasing ^2, 
can be understood as a superconductor-insulator transi- 
tion of charge 2e chargon (or Cooper) pairs. Note that 
a direct transition from the conventional insulator X to 
SC* is not generically possible. 

Figure M is a schematic phase diagram exhibiting the 
four phases - SC,SC*,T and T* - as well as the inter- 
vening transitions. Of the four, it is only in the band 
insulator X that spinous arc confined. In the other three 
phases the Z2 vortex is gapped out and uncondensed. 
These three phases exhibit excitations with "fractional- 
ized" quantum numbers. It is the condensation of the 
Z2 vortex which leads to confinement, leaving only the 
electron in the spectrum. 



1. Odd number of electrons per unit cell 



We now briefly consider the superconducting phases 
with odd integer filling, but still presuming s— wave pair- 
ing. Since chargon pairs are condensed in both SC and 
SC* , it suffices again to consider very large pair hopping 
amplitude, ^2- Moreover, with condensed chargon-pairs, 
the chargon operator can be replaced by the Ising spin, 
bj = Si = ±1 - the "ison" - as discussed above. After 
integrating out the gapped spinous, the effective theory 
again reduces to the Ising matter-plus-gauge theory as 
in Eqn. |l3^ , but with the addition of the Berry's phase 
term, Sb, 
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Note that the SC* phase is realized only for large K, as 
discussed above. In this limit, as we have emphasized 
several times, the effects of the Berry phase term Sb are 
expected to be innoccuous. Thus, SC* will continue to 
exist even in the presence of Sb- To see this in more 
detail, it is once again illuminating to pass to a dual rep- 
resentation, which exchanges the isons for the visons: 



FIG. 7. Schematic zero temperature phase diagram dis- 
playing the four phases SC, SC* , T, and I* . 
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The transition from SC* to SC occurs on condensing 
the ison - so that single chargons are themselves con- 
densed. Note that ison condensation leads to confine- 
ment of the excitations it has long-ranged statistical in- 
teractions with - the hc/2e vorton and the vison [i.e the 



with tanh{td) — e " and tanh{K, 
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Here ^ij 



is a dynamical Z2 gauge field, and as before ^1^ is an 
"applied" field with Y\^ ii^^* = — 1 through all spatial 
plaquettes. This theory is a direct Z2 analog of a U{1) 
superconductor in the presence of an applied magnetic 
field. 
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2. d-wave pairing and doping 

The discussion above generalizes readily to the case of 
d-wave pairing. In particular, a dSC* phase where char- 
gon pairs, but not single chargons, have condensed is an 
allowed phase in the model. It's properties are the same 
as for the s-wave case above, except that the spinous have 
a gapless d-wave dispersion. Also possible is a dSC phase 
coexisting with spin-Peicrls order, just as in the s-wave 
case. 

In the presence of finite doping with non-integer No, in 
either the s-wave or the d-wave case, the SC* phase is ex- 
pected to survive, since the Sb term is innoccuous in this 
phase. The conventional superconducting phases will be 
more sensitive to the value of Nq - several additional su- 
perconducting phases with broken lattice symmetries are 
presumably possible. 



FIG. 8. Schematic zero temperature phase diagram for the 
superconducting phases with an odd number of electrons per 
unit cell. The SC — SV phase is discussed in the text. The 
precise topology of the phase diagram when the couplings f c 
and K are both of order one is not firmly established. 



Consider briefly the phase diagram in the tc — K plane. 
A schematic phase diagram is shown in Fig. ^. Progress 
can be made in various limiting regimes. For Kd = 
the theory reduces to a pure Zi gauge theory with gauge 
field, liij = nl^^Hij. Since /ufj* plays no role in this hmit, 
the resulting phases are identical to that with even inte- 
ger A^o analyzed in the previous subsection. In particular, 
for large ic, we have a conventional superconductor SC 
with broken Z2 gauge symmetry, while for small ic, we 
get the exotic superconductor SC* . These phases survive 
for small Kd- It is easy to establish the absence of phase 
transitions for tc = 00 and AT = 0. For t^ = c», on the 
other hand, one can set /iy — fJ-ff^, and the model re- 
duces to the fully frustrated Ising model. As discussed 
extensively in Section IV, the results of Ref. p6t| show the 



existence of an ordered phase for large Kd where transla- 
tionally symmetry is spontaneously broken. In general, 
this is expected to lead to spin-Peierls order. In this case, 
though, the spin-Peierls order co-exists with supercon- 
ductivity. We will denote this phase as SC — SV. Sev- 
eral other ordered phases are presumably also possible 
though we will not discuss these here. 



In the SC — SV phase the external gauge field "pene- 
trates" with /iij Ri fJ-ij*, and the Ising model is frustrated. 
But as td is reduced, it eventually becomes favorable to 
"screen" out this external field, and enter a "Meissner" 
phase with (Jin My) ~ 1- When this happens the bro- 
ken translational symmetry disappears - along with the 
frustration - and one enters into SC. 



B. Vortex description 

In this subsection we show how the superconductor 
SC* may be described in the dual vortex formulation. 
The discussion in Section M was based on the action in 
Eqns. 109 for the spinous and hc/2e vortices. The sym- 
metries of the action allow the addition of "kinetic" terms 
for both Z2 gauge fields a and /i. Once again, though not 
present in the original action, these terms will be gener- 
ated upon integrating out high energy modes: 






(142) 
(143) 



It is of interest to explore the phase diagram for arbitrary 
positive values of the couplings Ka- and AT^ . We will show 
that the superconductor SC* emerges quite naturally for 
large K^r and A"^. As shown below, an important physical 
consequence of the addition of these Ka- and K^ terms is 
that the Chern-Simons term Scs is no longer effective in 
attaching flux to the vortices and the spinous. Note that, 
in the absence of flux attachment, the field e*^» creates 
a "naked" hc/2e vortex, i.e, an hc/2e vorton. Attaching 
flux of the field a, i.e a vison, converts this into a regular 
hc/2e BCS vortex. 

For ease of presentation, we specialize to the case of 
s-wave pairing and an even number of electrons per unit 
cell. In that case, the term Sb may be dropped from 
the action. Further, the spinous are gapped and can be 
integrated out. This will lead to an innoccuous renor- 
malization of the value of K^. 

In the vortex description, superconducting phases cor- 
respond to vortex vacuua. To analyse these, it is then 
appropriate to imagine integrating out the vortices. This 
will renormalize the value of K^ (or generate it if not 
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present originally). The resulting action has only the 
terms, 



S = Sa + Sa + Sn + ScS- 



(144) 



The term Sa leads to a gapless linear dispersing exci- 
tation (in the absence of long-ranged Coulomb interac- 
tions), and corresponds physically to the sound modes 
of the superconductor. The remaining three terms only 
involve the two Z2 gauge fields a and /i. As shown in Ap- 
pendix O, this action is equivalent to that of the Z2 gauge 
theory with Ising matter fields. If we choose to integrate 
out the /i, this is exactly the same as the Ising effective 
action derived in the previous subsection (Section VII A) 
to discuss the superconducting phases. Alternatively, we 
can inte grat e out the a field to obtain the dual theory as 
in Eqn. 14l|: 



O — Oil 



Su 



Jvis — J^a / J VjfJ-ijVj- 



(145) 
(146) 



Here tanh{K^) ~ e^^^" , so that K^ is the coupling dual 
to K^. Once again, Vi creates a vison, whose Z2 current 
is equal to the flux of the <t field. On the other hand, the 
vortex configurations of the gauge field n correspond to 
the ison excitations. 

As discussed earlier, the Z2 gauge theory with matter 
fields has two phases - a Higgs-confined phase and a de- 
confined phase. The Higgs-confined phase describes the 
conventional superconductor SC, and is perhaps easiest 
to understand in the limit in which both K^^ and K^ are 
small. With small X^, the gauge field is in its confining 
phase, so that test charges coupling to the gauge field /i 
are confined. There are actually two different particles 
minimally coupled to /i - the ^ vorton and the vison, 
with creation operator e*^* and Vi, respectively. As be- 
fore, the confined vorton-vison bound state is the con- 
ventional ^ BCS vortex. 

2e 

The deconfined phase describes the exotic supercon- 
ductor SC* . In this phase, test charges that couple to fi 
are deconfined. This implies that the || vorton and the 
vison are not bound together, and can propagate as inde- 
pendent gapped excitations, in agreement with the ear- 
lier discussion. In effect, within SC* the Chern-Simons 
term has been rendered ineffective and does not attach 
fiux. Also, configurations with vr flux in the gauge field /i, 
corresponding to the "ison" , exist as finite energy exci- 
tations. Thus, as before we conclude that there are four 
gapped excitations in SC* - the hc/2e vorton, the spinon, 
the vison and the ison. 

Note that a transition from SC* to an insulator ob- 
tained by condensing the hc/2e vortons (which are the 
fundamental U{1) vortices in this phase) leads naturally 
to the fractionalized insulator T* . This is because the 
vison is unbound from the hc/2e vorton in SC* , so that 



condensation of the latter leaves the former uncondensed. 
Indeed, the distinct excitations in the resulting insulator 
are the chargons, the spinous, and the visons - as appro- 
priate to X* . Thus, the exotic insulator I* may either 
he reached from SC by condensing — vortices or from 
SC* by condensing ^ vortons. In either case, the vison 
remains uncondensed. 

This completes the dual description of SC* . Compli- 
cations such as d-wave pairing or arbitrary filling Nq can 
be handled straightforwardly in this dual formulation as 
well, though we shall not do so here. 



VIII. EXTENSION AND GENERALIZATIONS 

A. General spatial dimension 

The Z2 gauge theory formulation (in the particle rep- 
resentation) is readily generalized to arbitrary spatial di- 
mension. The cases of physical interest are three dimen- 
sions (3d) and one dimension (Id), which we discuss in 
turn. For simplicity, we will restrict our attention to 
situations with integer filling per unit cell. The most 
important effect of spatial dimensionality enters in the 
properties of the pure Z2 gauge theory with action. 



S — Sa + Sb , 



(147) 



with Sb included when there are an odd number of elec- 
trons per unit cell. 



1. d^3 

In 3d and in the absence of Sb, the Z2 gauge theory 
again has two phases distinguished by the behaviour of 
the Wilson loop correlator ("area law" versus "perime- 
ter law"). As in 2d, the presence of Sb will enhance the 
stability of the perimeter phase, but the area law phase 
will still be present. The presence of the perimeter law 
phase, implies the existence of 3d insulators with electron 
fractionalization. But in constrast to 2d, the ffux tubes 
in the Z2 gauge field - the visons- are not point-like ex- 
citations, but become extended string-like excitations in 
3d. The area law phase again describes various confined 
insulating phases. Whether the presence of Sb leads to 
broken translational symmetry as in 2d is an interesting 
open question. Note, however, that in 3d it is not pos- 
sible to pass to a dual global Ising model. In fact, the 
pure Z2 gauge theory (in the absence of Sb) is in fact 
self-dual |13] in three spatial dimensions. 

To discuss the superconducting phases SC and SC* , it 
is necessary to understand the properties of the Z2 gauge 
theory coupled to Ising matter fields. In the absence of 
Sb, it is known pO| that in three spatial dimensions, 
there are again two phases - the Higgs-confined phase. 
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and the deconfined phase. These correspond to SC and 
SC* , respectively. Their distinguishing properties wiU be 
quaUtatively similar to the 2d case. As in 2d, we expect 
that the main effect of Sb would only be to make possible 
the existence of an SC phase with broken translational 
symmetry. 

In layered quasi- two dimensional systems, fractional- 
ized insulating phases in which each layer is decoupled 
from the others are possible, and exist as distinct phases 
from the isotropic ones discussed above. Such phases are 
currently under further investigation. 

Finally, it is worth emphasizing that while the exten- 
sion to 3d is straightforward in the particle representa- 
tion, the dual vortex representation necessarily involves 
string-like vortex degrees of freedom. 



2. d 



In one spatial dimension (Id), the Z2 gauge theory is 
always in its area law phase, with or without the Sb term. 
Thus, our formulation is incapable of describing electron 
fractionalization in one dimension. Evidently, fractional- 
ization in d = 1 must have different physical origins than 
for d > 1. To highlight this point, note that Id fractional- 
ization can be continuous, as exemplified by the spinless 
Luttinger liquid which supports charge-carrying excita- 
tions with essentially arbitrary (even irrational) charge. 
For d > 1, on the other hand, fractionalization is discrete 
- the fractionally charged excitations carry a definite ra- 
tional fraction of the electron charge. As in the fractional 
quantum Hall effect, this discreteness can be traced to the 
binding (and condensation) of a discrete number of vor- 
tices. This physics appears to be qualitatively different 
than the "solitonic" mechanism responsible for fraction- 
alization in Id, 



B. Finite temperature 

In our formulation there is a sharp distinction between 
fractionalized and confined phases at zero temperature, 
which is independent of whether or not the phases in 
question have any sort of conventional long-ranged order. 
It is extremely interesting to ask whether this sharp dis- 
tiction survives at finite temperature. Consider first the 
deconfined phases in 2d. In these phases, the point-like 
vison excitations are gapped at zero temperature. How- 
ever, since the energy cost to create a vison is finite, at 
any non-zero temperature there will be a non-vanishing 
density of thermally excited visons. In the absence of 
other kinds of order (eg. magnetic), this low tempera- 
ture regime will be smoothly connected to the high tem- 
perature limit, without an intervening finite temperature 
transition. Thus, in 2d the sharp distinction between 



fractionalized and confined insulators does not survive at 
finite temperature. 



But in 3d, the vison excitations in the deconfined phase 
are string-like extended objects, with an energy cost pro- 
portional to their length. Consequently, at low temper- 
atures arbitarily large vison loops will not be thermally 
excited - the vison loops will be "bound". As the tem- 
perature increases, there will be a transition at which 
the vison loops unbind and proliferate. Thus, the frac- 
tionalized insulator in three spatial dimensions undergoes 
a finite temperature phase transition associated with the 
unbinding of vison loops. A defining characteristic of the 
low temperature phase is that vison loops will cost a free 
energy linear in their length. Equivalently, hc/2e (or Z2) 
magnetic monopole "test charges" are confined even at 
finite temperature, with an infinite free energy cost to 
separate them. A confinement of monopoles is also one 
of the characteristics of a 3d superconductor, but quite 
remarkably the confinement here is occuring in a "nor- 
mal" non-superconducting phase. The conventional in- 
sulating phases with confinement at zero temperature, 
on the other hand, will not exhibit finite temperature 
transitions (other than those associated with the loss of 
conventional long-ranged order - eg. magnetic). 



To understand the origin of these results, we briefly 
discuss the properties of the pure Z2 gauge theory (with 
no matter fields) in 3 -I- 1 space-time dimensions in more 
detail. At zero temperature the theory is self-dual |lj] 
- the duality transformation interchanges the "electric" 
and "magnetic" fields of the gauge theory. For K > Kc 
when the gauge theory is in its deconfining phase, the 
theory has string-like vison excitations (which are Z2 
"magnetic" flux tubes) with a finite energy cost per 
unit length. For K < Kc the gauge theory confines 
with area law Wilson loops, but there are nevertheless 
string-like excitations in this phase as well. These can 
be understood via duality, which interchanges the area 
and perimeter law phases - the string-like excitations in 
the area law phase are simply flux tubes of the dual Z2 
gauge field. Physically, these dual tubes are "electric 
flux tubes" responsible for the confinement of "electric" 
charge in the area law phase. Specifically, when two test 
Z2 "electric" charges separated by a distance R are in- 
troduced into the system, the resulting "electric" flux is 
concentrated in a tube that extends from one test charge 
to the other with an energy cost proportional to i? - the 
linear confinement. Similarly, in the perimeter phase, 
dual test charges {Z2 "monopoles") that act as sources 
for the visons are confined. 
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FIG. 9. Schematic finite temperature phase diagram for 
the pure Z2 gauge theory in three spatial dimensions. Upon 
including the coupling to the chargons and the spinous, the 
finite temperature transition for K < Kc is smeared and be- 
comes a crossover only, while the one for K > Kc continues 
to exist. 



Now consider the properties of the gauge theory at fi- 
nite temperature. The phase diagram is well-known pM 
and is shown in Fig. BI There are three finite temper- 
ature phases. For K > Kc, at small but non-zero tem- 
peratures, large ("magnetic") vison loops are bound as 
their energy cost is proportional to their length. Simi- 
larly, for K < Kc at low temperatures, large "electric" 
flux loops are bound. At high temperature, for any K , 
both kinds of loops are unbound. The transition from the 
low temperature to the high temperature phase is there- 
fore associated with the unbinding of (electric) magnetic 
vison loops for K (lesser) greater than Kc- 

In the low temperature phase for K < Kc, the free en- 
ergy of an isolated static test "electric" charge diverges, 
so that test charges are confined. In the high temperature 
deconfined phase, the free energy cost is finite. Formally, 
the pure Z2 gauge theory has a global Ising symmetry at 
finite temperature which is broken in the high temper- 
ature phase. As shown by Polyakov [pl| , a convenient 
characterization of this transition is through the follow- 
ing operator: 



Lr 



M-1 
n=0 



"^r-^-iiT ,r-\-(n-i-l)T 



(148) 



for each site f of the spatial lattice. Here f is a vector 
along the (imaginary) time direction of length the time- 
slice. The product is over all the temporal links at that 
site, and M is the number of time slices. This operator 
Lr is often referred to as the "Polyakov loop" . The free 
energy F{r,r') to introduce two test charges at sites r, 
r' is directly related Bit] to the correlator of Lr through 



-F{r,r')/T _ 



<^ L/r^r' -^ 



(149) 



Thus, test charges will be confined if this correlator goes 
to zero at large distances - on the other hand, if this 
correlator goes to a constant, the test charges will be 
deconfined. Furthermore, consider the following trans- 
formation on the gauge fields 



ap. 



r+nof ,r+(no + l)T 



eap. 



r+noT,r+{no + l)7 



(150) 



where e = ±1 independent of r, and uq is fixed. The 
action of the pure gauge theory is invariant under this 
transformation - implying a global Ising symmetry of the 
theory. The operator Lr, however, transforms as 



eLr 



(151) 



Thus Lr is an order parameter for this global Ising sym- 
metry. In the low temperature phase for K < Kc, Lr has 
no expectation value, the global Ising symmetry is unbro- 
ken, and test charges are confined. At high temperatures, 
however, Lr acquires an expectation value breaking the 
global Ising symmetry, and the test charges are decon- 
fined. 

For K > Kc, the self-duality of the Z2 gauge the- 
ory implies the existence of a dual global Ising symme- 
try, with an order parameter that is the dual analog of 
the Polyakov loop. In the low temperature phase, this 
dual global symetry is unbroken - in this phase dual test 
charges (i.e Z2 monopoles) are confined. At high tem- 
peratures this dual global symmetry is broken and the 
dual test charges are deconfined. 

Consider next the effects of coupling matter fields (the 
chargons and the spinous) to the Z2 gauge field. As these 
carry Z2 gauge "electric" charge, it is easy to see that the 
action is n o longer invariant under the transformation in 
Eqn. 
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Indeed, this transformation is equivalent to 
changing the boundary conditions on the chargon fields 
from (/?— )periodic to anti-periodic, and vice versa for 
the spinous. Moreover, if the matter coupling is weak, 
the matter fields may formally be integrated out |^^ to 
leave behind a "magnetic field" term that couples linearly 
to the Polyakov loop order parameter of the global Ising 
symmetry. There is then no longer any transition sepa- 
rating the low and high temperature regimes. Physically, 
this is exactly as expected - for i^ < Kc, the electronic 
system is in a conventional confined insulating phase at 
zero temperature. 

On the other hand, since the chargons and spinous do 
not carry any dual Z2 "magnetic" charge, the dual global 
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Ising symmetry rem.ains even in their presence. The finite 
temperature transition for K > Kc should thus remain in 
tact. Consequently, we arrive at the striking conclusion 
that the three dimensional fractionalized insulator un- 
dergoes a finite temperature transition associated with 
the unbinding of vison loops. This conclusion will not 
be affected by the Berry's phase term Sbi which is quite 
innocuous in the fractionalized insulator. 



C. Spin-rotation non-invariant systems 

The Z2 gauge theory formulation (in either the parti- 
cle or vortex representations) works equally well in the 
absence of spin rotation invariance. In particular, frac- 
tionalized phases continue to exist even when spin is not 
a good quantum number. (Spinless fermion systems can 
also be handled with no fundamental modifications). For 
these reasons, we have avoided the term "spin-charge sep- 
aration" , in favour of the more general term "electron 
fractionalization" . 



D. Analogies with nematics 

Certain aspects of our formulation might be familiar 
from the classical statistical mechanics of nematics. The 
order parameter for a nematic is a headless three com- 
ponent vector. Lattice models of nematics are usually 
formulated in terms of an ordinary three component vec- 
tor - the headless nature being incorporated through a 
local Z2 gauge symmetry which inverts the local vector 
order parameter. Here, we briefly explore the analogies 
between the classical phases of nematic systems, and the 
quantum phases discussed in this paper. 

The analogy is closest if we consider s-wave pairing 
with an even number of electrons per unit cell, and fur- 
ther, integrate out the spinous to work with just the char- 
gons and the a field. The action describing the system is 
then, 



S = — 2ic / cTiji 



;)-^En^. 



(152) 
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As formulated, this describes a gitanium problem of char- 
gons coupled to a fluctuating Z2 gauge field in two spatial 
dimensions. But alternately, we may view it as a classical 
Hamiltonian for a three dimensional XY nematic. In- 
deed, an 0(3) version of the same model was introduced 
a few years ago by Lammert, Rokhsar, and Toner pSl to 
describe nematic ordering in three dimensions. Further, 
they argued that their lattice gauge nematic model ad- 
mits three distinct phases - an ordered nematic phase, 
and two isotropic phases. The nematic phase breaks the 
rotational symmetry, and the Z2 gauge symmetry. For 
an XY system, this is the direct analog of the supercon- 
ducting phase. Moreover, the physical hc/2e vortices of 



the superconductor correspond directly to the "disclina- 
tions" in the nematic fluid. 

The two isotropic phases in the nematic are distin- 
guished Q by the free energy cost per unit length to 
externally impose a disclination line through the system. 
In particular, in the conventional isotropic phase, the free 
energy cost per unit length is zero (as the length goes to 
infinity). The disclinations are condensed. But, in the 
unconventional isotropic phase 133], the free energy cost 
per unit length is a constant (as the length goes to in- 
finity). In the context of this paper, the isotropic phases 
correspond to insulating phases. As we have elaborated 
at length, there are two insulating phases X and I* which 
are distinguished by whether or not the visons (which are 
the relics of the hc/2e vortices in the insulating phases) 
are condensed. Thus, the conventional insulator cor- 
responds, in the nematic analogy, to the conventional 
isotropic phase. Note that the energy cost of a vison 
(which is the action cost per unit length of the world- 
line) is zero in this phase. Similarly, the fractionalized 
insulator 2* corresponds to the unconventional isotropic 
phase of the XY nematic. In I* the visons have finite 
energy cost, again just like the disclination lines in the 
unconventional isotropic fluid. 

The phase transition between SC and either insulating 
phase is second order. In contrast, for the 0(3) nematic 
system considered in Ref. p3 , the transition between the 
nematically ordered phase and the conventional isotropic 
phase is first order, while that to the other isotropic phase 
is second order. This difference is due to the XY sym- 
metry of the superconducting system, as opposed to the 
0(3) symmetry of the nematic. 

For the more general situation, with coupling to the 
spinous or with an odd number of electrons per unit 
cell, a direct correspondence with the nematic system 
no longer holds. Nevertheless, we believe that the dis- 
cussion in this subsection may help (some) readers get 
further intuition and insight into our formulation. 



IX. RELATION TO PREVIOUS APPROACHES 

We now comment on the connection between the Z2 
gauge theory and earlier approaches to electron fraction- 
alization. We begin by making contact with earlier pa- 
pers on the "nodal liquid" . Earlier formulations of the 
nodal liquid (in Ref. M and [EOJ) focusscd on the impor- 
tance of "vortex-pairing" as a means to describe charge 
fractionalization in two-dimensions. In Ref. M a the- 
ory was formulated in terms of vortices in a local super- 
conducting pair field, and shares many features with the 
approach taken here, particularly the dual formulation 
detailed in Section M. In Ref. flG], Chern-Simons theory 
was used to convert spinful electrons into bosons, and a 
dual formulation was developed in terms of vortices in 
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these bosonic fields. The Z^ gauge theory and its dual 
Ising Chern-Simons vortex theory developed in this pa- 
per, not only ties together both earlier approaches into a 
unified framework, but allows for a more direct quantita- 
tive analysis of "microscopic" models. We now describe 
this connection in a bit more detail. 

In Ref. PI a spinon operator was defined as an electron 
with its charge screened by "one-half of a Cooper pair. 
The latter coincides precisely with the "chargon" intro- 
duced in Eqn.E9^, showing the equivalence of the spinous 
as well. The importance of the long-ranged interaction 
between the spinon and hc/2e vortex was emphasized in 
Ref. 0. It was suggested that this interaction could be 
implemented by employing a U{\) Chern-Simons term 
to attach flux to both species of particles. But since the 
spinon number is not conserved, it was suggested that the 
flux could be attached to the (conserved) z-component of 
the spin. Moreover, it was argued in Ref. ^ that due to 
the statistical interactions, condensation of hc/2e vor- 
tices should lead to confinement of spinous. In the dual 
vortex formulation presented in this paper the statisti- 
cal interaction between vortex and spinon is described 
in terms of a novel Ising-like Chern-Simons term. It is 
important to stress that this does not require the spin of 
the spinon to be conserved, in contrast to the U{\) ap- 
proach, since the "Ising-flux" is attached to the conserved 
Z2 charge of the spinous. Moreover, the Ising formulation 
clearly shows that condensation of the hc/2e vortices - or 
the visons - leads to confinement of spinons and chargons. 
In the global Ising model for the visons with {vi) ^ 0, the 
linear confinement is due to the required line of negative 
Ising couplings connecting the two spinons. In the Z2 
gauge theory formulation, it follows from the area law 
for the Wilson loop. 

In Ref. |lO|| , a theory was developed by converting spin- 
ful electrons into spinful bosons - using Chern-Simons to 
attach flux to the electrons spin - and then passing to 
a dual representation of vortices in these bosonic fields, 
denoted $q, with spin label a =t, j. A lattice version 
of this theory can be written in terms of the phases, 9 on 
of the vortex field operators, $„ = e*^", with effective 
Euclidian action. 
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cos{6i, 



+ a") + 5c.K). (153) 



Here, ?, j label sites of the 2-1-1 space-time lattice, t^ is a 
dimensionless vortex "hopping" term and Scs is a Chern- 



Simons terms involving the field a" 



iK The curl of 



a" corresponds to the conserved electrical current of the 
electrons with spin a. In Ref. [^ , two different composite 
"pair" vortex operators were considered; 



$p = $T$i = e'" 



$^ = $t$| =e'' 



(154) 



which are minimally coupled to af'"^ = a^ ± o^ , respec- 
tively. The action can be re-expressed in terms of these 
composite phase fields using the relation 
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Here, the Ising spins Vi = ±1 are the "visons". The pri- 
mary emphasis of Ref. |1^ was an analysis of fractional- 
ized phases, such as the nodal liquid. It was emphasized 
that fractionalization occurs when (vi) = 0, and breaking 
the Ising symmetry with (vi) ^ corresponds to confine- 
ment. Deep within the deconfined phase it is possible to 
integrate out the massive visons, which generates local 
terms such as, 
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which 
and 



SkcI. = -t2.cos{Ql - 0; + a%), (157) 

ich describes the hopping of the hcj e vortex pair, $p, 
1 



S. 



spinon 



-t.cosiQ, — 
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Due to the Chern-Simons terms above, this corresponds 
to the hopping of fermionic spinons which carry Sz = 1/2. 
The relationship between this formulation, in terms of 
"electron" vortices, and the dual vortex theory of Section 
constructed in terms of BCS hc/2e vortices is at first 
not apparent. But consider introducing a vortex opera- 
tor, $ — e*^, whose square equals the hc/e vortex pair 
operator: $^ — $„. This requires that. 



\Op + ^{l-v), 
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which implies. 



$ = i-e'^-/^ 



(159) 



(160) 



As defined $ carries vorticity hc/2e, and can tentatively 
be identified as the BCS vortex. To complete this identi- 
fication it is necessary to show that there is a long ranged 
statistical interaction between this hc/2e vortex and the 
spinon. Evidence for this is provided by the following 
argument. We first imagine explicitly adding the vortex 
hopping term Sspinon to the action in Eqn. HSfl. We then 
absorb the field 0°" into a" . We may now re-express the 
action Eqn. 



156 in terms of 9j 



S 



with 



-ty'^^ijCos{9i - 9j 



1 

— a, 
2 ■ 



IJ-v 



cos{-a'^j). 



(161) 



(162) 
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Here, we have defined a^ 



vortex hopping term Sspinon above, if we specialize to the 
hmit of large t^, it is legitimate to restrict of- to be 27r 
times an integer. With that restriction the gauge field 
/i.y = ±1, reducing to an Ising Z2 gauge field. Now, 
imagine putting a stationary spinon on one site of the 
original spatial lattice. In this dual vortex representa- 
tion this corresponds to a plaquette with A x a'^ = 27r, 
or equivalently to a product fl^ ^ij = — 1 for all plaque- 
ttes pierced by the spinon "world-line" . Since the hc/2e 
vortex is minimally coupled to Hij, this establishes that 
it does indeed acquire a minus sign upon being trans- 
ported around a spinon. In the dual vortex formulation 
in Section |V|, a tt-Aux tube in fiij is attached to each 
spinon by the Ising-like Chern-Simons term. To complete 
the mapping between these two formulations requires, fi- 
nally, to re-fermionize the spinon creation operator, e*^» , 
(fermionic due to the Chern-Simons term 5*28 [a"^]) effec- 
tively replacing it with spinful fermions fia- 

Finally we comment briefly on the relationship with 
theories based on slave boson/fermion approaches to elec- 
tron fractionalization. A number of authors have exam- 
ined insulating Heisenberg antiferromagnetic spin mod- 
els in the hope of finding phases with deconfined spinon 
excitations through these approaches. However this pro- 
gram has generally been quite unsuccesful - the U{1) or 
SU{2) gauge symmetry introduced in the slave boson or 
fermion representations ultimately leads only to confined 
phases. A notable exception however is the work of Read 
and Sachdev |13| on large- A^ Sp{2N) frustrated antiferro- 
magnets, and related quantum dimer models |17|. Under 
certain special conditions, these authors demonstrated 
the existence of quantum disordered phases with decon- 
fined spinous in their theory. It is worth pointing out that 
fractionalization is achieved when the U{1) gauge sym- 
metry (introduced by the Schwinger boson representation 
of the Sp{2N) spins) is broken down to Z2 by conden- 
sation of pairs of bosons. The fully frustrated transverse 
field Ising model appears in that description as well ^^ . 

Slave boson representations of electron operators have 
been used extensively to discuss spin-charge separation 
issues in doped t — J models. However, the resultant 
compact [/(I) or SU{2) gauge theories presumably al- 
ways lead to confinement, unless the gauge symmetry is 
broken down to Z2. This may be be achieved by pairing 
the spinons |I^. Indeed, the slave-boson mean field treat- 
ments of the t— J model do find pairing of spinons below 
a finite temperature at low doping. As we have empha- 
sized in this paper though, even in the undoped limit and 
without frustration, the Heisenberg spin model may be 
rewritten in terms of fermionic spinon operators coupled 
to a fluctuating Z2 gauge field. Equivalently spinon pair- 
ing terms may be added to the Hamiltonian describing 
the Heisenberg magnet without altering any of the phys- 
ical symmetries. We have shown that electron fractional- 
ization is definitely possible once charge fluctuations are 



In the presence of the incorporated into the description. 



X. CONCLUSION AND DISCUSSION 
A. Summary 

The primary focus of this paper was to explore the pos- 
sibility of electron fractionalization in strongly correlated 
electron systems in spatial dimension greater than one, 
and in the presence of time reversal symmetry. We based 
our discussion on a particular class of microscopic models 
designed to capture the physics essential to the cuprates, 
although our decription of fractionalization is more gen- 
eral. Starting from these models, we developed a new 
gauge theory of strongly correlated systems consisting of 
charge e, spin-zero bosons (the "chargons") and charge 
zero, spin 1/2 fermions (the "spinons"), both minimally 
coupled to a fluctuating Z2 gauge field. Remarkably, the 
spin-sector of the theory at half-filling and in the absence 
of charge fluctuations, is formally identical to a spin one- 
half Heisenberg antiferromagnet. In this limit the Z2 
gauge field enforces the constraint that the spinon num- 
ber on each site is odd - physically equivalent to the single 
occupancy constraint, imposed with additional unneeded 
redundancy in earlier U{1) gauge theory formulations of 
the Heisenberg model. 

Charge fluctuations, however, are naturally incorpo- 
rated into our Z2 gauge theory, and when they be- 
come large the theory describes a dxi_yi superconduc- 
tor. Analysis of the theory in the intermediate region 
reveals that there are two qualitatively different routes 
for the evolution from the antiferromagnet to the super- 
conductor. One route is through conventional insulating 
phases in which fluctuations of the Z2 gauge field con- 
fines together the chargon and the spinon, leaving only 
the electron in the spectrum. But a more interesting 
possibility takes one through phases in which the elec- 
tron is fractionalized, and the chargons and spinons exist 
as deconfined excitations. With dx2_^y2 pairing, this frac- 
tionalized insulator is the nodal liquid [gl|,^ , with gapless 
spinon excitations at four points of the Brillouin zone. It 
seems likely that the ultimate transition from the insu- 
lating phases to the d^^-y-^ superconductor occurs close 
to the boundary between the confined and deconfined in- 
sulating phases. Thus, which of these two qualitatively 
different routes is realized in any particular experimental 
system could depend sensitively on microscopic details. 

In addition to the chargons and spinons, the 2d nodal 
liquid supports Ising-like point excitations - the "visons"- 
which correspond to vortices in the Z2 gauge field. These 
gapped vison excitations play a central role in our anal- 
ysis of fractionalization, as becomes clear upon passing 
to a dual description in terms of hc/2e BCS vortices (of 
a conventional superconductor) and the spinons. In this 
dual framework, the nodal liquid can be accessed by a 
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pairing and condensation of the hcjle vortices, as em- 
phasized in ealier work |^Jl^ . This reveals that the vison 
excitations are simply the remnant of the unpaired hc/2e 
vortices which survive in the insulating nodal liquid. 

The utility of the vison excitations goes far beyond 
giving a simple description of the nodal liquid. Indeed, 
the pure Z2 gauge theory in 2 + 1 space-time dimensions 
is dual to the global 2-1-1 dimensional Ising model - 
and the Ising spins are simply the vison creation oper- 
ators. Remarkably, an unusual Berry's phase term in 
the gauge theory corresponds simply to frustration in 
the dual Ising model, with full frustration at half-filling. 
The fully frustrated quantum Ising model arose in earlier 
work by Sachdev and coworkers |ll6|,y_7| in their analysis 
of frustrated magnets. Ordering the dual Ising model by 
condensation of the visons, generally will break transla- 
tional symmetry and lead to conventional confined in- 
sulating phases such as the spin-Peierls phase. In three 
spatial dimensions (3d) , the visons become loop- like exci- 
tations, and are closely related to vortex-line excitations 
which occur in a conventional superconductor. Surpris- 
ingly, this implies that a 3d fractionalized insulator "sur- 
vives" at finite temperature, being separated from the 
high temperature regime by a finite temperature phase 
transition. As in a conventional superconductor, the 3d 
fractionalized insulator confines hc/2e monopole excita- 
tions even at non-zero temperature. 

Within the Z2 gauge theory approach, a conventional 
superconductor is described as a condensate of charge e 
chargons. A superconducting phase involving condensa- 
tion of chargon pairs (z.e Cooper pairs) without conden- 
sation of single chargons was shown to exist - this has 
several exotic properties distinguishing it from the con- 
ventional superconductor. 



B. Experiments 

We close with a very brief discussion of some of the 
experimental signatures of electron fractionalization. As 
we will see, experimental detection of fractionalization 
may be quite subtle. Further theoretical understanding 
of fractionalized phases leading to detailed experimen- 
tal predictions are clearly called for. Our discussion will 
necessarily be brief. 



1. Two dimensional nodal liquid 



Earlier work on the nodal liquid |21^,P[ outlined a num- 
ber of experimental signatures of the two-dimensional 
nodal liquid, and we have little to add here. As pointed 
out in the earliler papers, perhaps the most telling indi- 
cation will be in angle resolved photoemission (ARPES) 
which directly measures the electron spectral function as 
a function of the momentum fe, and frequency to. As the 



electron is fractionalized into the chargon and the spinon 
in the nodal liquid, its spectral function will not have a 
sharp quasiparticle peak even at zero temperature. Note 
that bound states of the chargon and the spinon (which 
could lead to sharp spectral features) are not expected 
here at low energies as the spinons are gapless. 



2. SC* 

We have discussed the basic physics of the exotic su- 
perconduc tor S C* obtained by condensing chargon pairs 
There are several qualitative experimen- 



in Section VII 



tal distinctions between this phase and the conventional 
superconductor which we now briefly discuss. The most 
striking is again in the electron spectral function as mea- 
sured in ARPES. As discussed in Section VII, the elec- 



tron decays into a spinon and an Ising part of the charge 
- the "ison" excitation. Thus, we expect that the elec- 
tron spectral function does not have a sharp quasiparticle 
peak in the SC* phase. Again, since the isons are massive 
excitations while the spinons are gapless, bound states 
of the two are generally not expected at low energies. 
The presence of gapped ison excitations would also affect 
the thermodynamics, and contribute to the thermal con- 
ductivity at some intermediate temperatures. However, 
these signatures are likely to be quite subtle. A striking 
theoretical feature of SC* is that the conventional BCS 
hc/2e vortices are splintered into pieces - the U{1) "vor- 
ton" carrying the circulating electrical currents, and the 
Z2 vison. Since the spinons do not have a long-ranged 
statistical interaction with the hc/2e vorton, it is tempt- 
ing to speculate that the structure of the core states in 
such a vorton would be qualitatively different from that 
of an hc/2e vortex in a conventional superconductor. 



3. Three dimensional effects 

In striking contrast to a two dimensional nodal liquid, a 
genuinely three dimensional nodal liquid has a finite tem- 
perature phase transition associated with the unbinding 
of vison loops. This phase transition could lead to ob- 
servable singularities in the measured properties of the 
system. But due to the highly anisotropic nature of the 
cuprates, it is perhaps more natural to speculate that a 
fractionalized phase would consist of decoupled 2d sys- 
tems, with a confinement of spinons within each layer. 
Clarification of such interlayer confinement physics will 
be necessary in order to disentngle the subtle interlayer 
behavior of the cuprate materials, both in the normal 
and superconducting phases. 
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APPENDIX A: PATH INTEGRAL 

In this Appendix we derive a path integral expression 
for the partition function of the spinon-chargon Hamil- 
tonian. A crucial role is played by the constraint on the 
Hilbert space, which naturally introduces an Z2 gauge 
field. 

To this end, we work with fermionic coherent states 
built from the spinon operators, /„ and /^, which are 
defined in the standard fashion: 



l/c) = e 



(/a 



- -Ufl 



|0), 



p/q/q 



(Al) 



(A2) 



where the spinon operators are denoted with "hats" , and 
fa and fa are Grassman numbers. The bra and ket states 
denoted with a "0" , are fermionic fock states with no 
spinous present. Here we have suppressed the depen- 
dence of the fermion operators and Grassman fields on 
the spatial coordinate, r. In the charge sector of the the- 
ory we choose a basis of states diagonal in the phase of 
the chargon field, denoted \(j)). 

The partition function in Eqn. |4^ can then be ex- 
pressed as 



Z ^ dfadfa 



2Tr 



-fcfc 



{-f^;^\(e-^"Vr\fa-A) 



(A3) 

with e = 0/M and V the projection operator defined in 
Eqn. p9l Inserting the resolution of the identity between 
each time slice gives, 



^ = n / dfradfrc.d^re~^'-^--'Mr, 



T = l 

with matrix elements 



Mr^{fr;(br\e~'''r\fr;(br-l) 



(A4) 



(A5) 



and appropriate boundary conditions on the fields: 
/m+1 = -/i and 4>o = 4>m- 

The matrix elements can be readily evaluated for small 
e by inserting a complete set of states diagonal in the 
chargon number, N . Using the definition of the projec- 
tion operator in Eqn. B3 gives, 



_. 00 

(A6) 



crT=±l N^ = — oo 



with 



Er = e 



_ -eH{N^,4,^.U,a^U) 



(A7) 



Upon making the change of variables in the Grassman 
functional integral, 

^r,/r ^ fr, (A8) 

the full partition function can finally be re-expressed as, 

/M 00 
l[dfrdfrdc^r ^ J^ S"^' ^^9) 

.. — 1 i\r — ^^ — LI 



N^ = -oo ar = ±l 



with, 



M 



S=Sl + St + eY,H{Nr,<jyrJrfr). (AlO) 



with 



M 



Sl = ^[frif^T+lfT+l - /r)], 



(All) 



T=l 



and 



M 
St^-iY,Nr[<Pr-<Pr^l + ^{l-CTr)]. (A12) 



T = l 



Throughout, we have suppressed the explicit r and a 
subscripts on the fields, displaying only the time-slice 
dependences. 



APPENDIX B: Z2 GAUGE THEORY WITH 
0^2 _Y-^ PAIRING 

In this Appendix, we will provide the outline of a mi- 
croscopic derivation of the Z2 gauge theory in the pres- 
ence of dr^2_y2 pairing correlations. We begin with the 
Hubbard-type Hamiltonian Eqn. discussed in Section 



H = Ho + Hj + Ha + Hu 



(Bl) 



The crucial difference with the s-wave case is in the struc- 
ture of the "pairing" term H/^. 

We now follow exactly the same strategy as in the s- 
wave case, defining chargon and spinon operators. A path 
integral representation of the partition function is readily 
set up with the main difference being in the pairing term 
which becomes 



32 



SA^e Yl ^rr'{Kbr'+C.c)Brr', (B2) 

B„' = Kr'Ur^Ir'i - (T^i) + C.c). (B3) 

We have suppressed the r index on all fields. It will be 
convenient to use a slightly different decoupling of the 
Hj term. We write 



APPENDIX C: ISING SELF-DUALITY 



In this Appendix, we will review the self-duality of the 
Z2 gauge theory with matter fields in 2 + 1 dimensions. 
As a limiting case, we recover the duality of the pure Z2 
gauge theory to the global Ising model. The theory is 
defined by the lattice action 



[dXrr'dXrr'dVrr'dr]*r,]e ^"^ 



(B4) 
(B5) 



Shs = Shs[x] + ShsM, 

Shs[x]^<^ X! J['^\Xrr'\'^ - iXrr'frafr'a+C.c)], (B6) 

Shs[v]^e J2 A'^-nw? (B7) 

Krr'y .r 

+ (rjrr'arr'ifrhfr'l - Irjr']) + C-c)- (B8) 

Here arr' = +1 for bonds along the x-direction, and 
equals —1 for bonds along the y-dircction. Note that 
'S'fts[x] is the same as before. This decoupling of the spin- 
spin interaction is standardly used in the SU{2) gauge 
theory formulations of the t — J model. We emphasize 
though that our formulation has, as we will show, only 
an Z2 gauge symmetry. We now shift the two Hubbard- 
Stratonovich terms: 



S[s,(t] = Ss + Sa, 



(CI) 

(C2) 
(C3) 



The constants J, K are assumed to be positive. The in- 
dices i, j label the sites of a three dimensional cubic lat- 
tice. It is convenient to first rewrite the SiGijSj term or 
each bond using the following identity: 



j-,'> Si (7ij Sj 



A ^ exp [2 J, 



d'^ij 



rii,-— 0,1 



• " / -1 N 



Xrr' ^ Xrr' j^r"r' ■, 

rjrr' -^ Vrr' + -jiK^r' + C.c) . 



Here tanh{Jd) 



= e-2^. 



and A = 



(B9) 
(BIO) 



V^- 



(C4) 
(C5) 
From now on, 



The shift of x is as before, and eliminates the spinon- 
chargon interaction coming from rewriting the electron 
hopping term. The shift of 77 eliminates the pairing term. 
The net spatial part of the action is then. 



Sr^eY, 2J{\Xrr'? + l^rrf ) + ^cr + S], + S^, (Bll) 
Krr'y 

Scr = -e J2 [(2tXrr' + 2A(77„, + V*rr'))Kbr' + C.c], (B12) 

^Ir ^ ^'^ Y JXrr'frafr'a + C.C, (B13) 



we will drop the constant A as it just contributes to an 
overall multiplicative constant to the partition function. 
The Uij take the values 0, 1. Upon using this identity for 
every bond, and doing the sum over Si, we get 

exp(-^,) = rr.,,Tr„.^. (j[cos (f (A.n)) j (C6) 

exp 2Jd Y ^^3 + Yl * 2"«j(l " '^^3) ■ C^^) 

\ ij i / 



Here A.n is the lattice divergence of the link variable n. 
We now notice that the cosine can be written as 



^sr — Vrr' ^rr' [fr^ fr' i — Irlfr']) + C.C. 



(B14) 
(B15) 



(A.n) =(-1)^5 (-l)^-M 



(C8) 



The shift in -q also generates a "Cooper pair" hopping 
term cos(2(j)r — 2(f>r') with a negative hopping amplitude 
of order b? j J . This is not expected to be important for 
the issues of fractionalization that we primarily wish to 
discuss. So we will for the most part drop it. 

The X, 77 integrals may be done by saddle point - a 
uniform, real saddle point solution < Xrr' >= Xvn < 
r\„i >— r]o breaks the Z2 gauge symmetry. Parametriz- 
ing the fluctuations about it by Xrr' — Xo<^ij , Vrr' — Vo<^ij 
as before, we arrive at the Ising gauge theory appropriate 
for the dx2_y2 superconductor. 



where S{m, n) is the Kronecker delta function for two in- 
tegers TO, n. The term multiplying the delta function is 
a total derivative that contributes zero on summing over 
all sites - we will therefore drop it. Note that the delta 
function imposes conservation modulo 2 of the link vari- 
able Uij at every site. This conservation can be made 
more explicit by defining a Z2 current a: 



(~ir 



(C9) 



We now solve the current conservation condition by writ- 
ing the Z2 current a on any link as the fiux of a dual Z2 
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gauge field /i through the plaquette of the dual lattice 
pierced by this link: 



a. 



(-!)"-■ =n^'^- 



(CIO) 



The fiij are understood to be defined on the links of the 
dual lattice, and the plaquette product for the /x is around 
the appropriate plaquette of the dual lattice. Note that 
this is directly analogous to the standard duality trans- 
formation of the XY model. 

We next solve for the Uij in terms of the /iy : 



1 - IId ^^^j 



(Cll) 



The Uij may now be eliminated from the action in favor 
of the fXij . The result (after dropping overall multiplica- 
tive constants) is the following identity 

^/E.-^'--^^^ ^^exp(-5^-5c5), (C12) 

Si fl 

S^^JdY^Hfi^j, (C13) 

<U> \ n / 

The last term has a structure similar to a Chern- 
Simons term, but for the group Z2. It's exponential is 
actually invariant under cr ^^ /i. This can be seen as 
follows. Write 



-Scs 



n nn 



(C15) 
(C16) 



In the last equation. Ax fi is the lattice curl of /i on the 
plaquette of the dual lattice pierced by < ij >. If we 
now perform a lattice integration by parts, we get 

= exp I - ^ ^^ y'^'') (1-^^^)) ' (^18) 

where now the sum is over links < ij > of the dual lattice. 
The full partition function can then be written as 



Z = Tr^^^ exp {-S^ - 5"^ - Scs) 



(C19) 



The duality of the full action is now apparent. In par- 
ticular, the action is invariant under the exchange cr <-> /u, 



Jd *^ K. To make the duality even more explicit, we 
again use the identity Eqn. C12 to write 



y^ exp {-Sa - Scs) = E °^P ^<^ E ^«My"j 

a Vi \ ij 

(C20) 

where v-i — ±1 and tanh{Kd) = e^^^ . The partition 
function now becomes 



Z 



= Trr e^''^'j'''''*'"'"^''''^°nn 



Mij 



(C21) 



which is exactly of the same form as in terms of the 
original variables {si,aij), but with the dual couplings 
{Jd, Kd), thus establishing the self-duality of the theory. 
As a special case, consider the limit when J = 0. Then 
the action in Eqn. CI is that of the pure Z2 gauge the- 
ory. Under the duality transformation, we now get the 
form Eqn. C21 but with the dual coupling Jd — cxd. This 
means that the fluctuations of the dual gauge field /i are 
frozen - we may choose a gauge in which fXij = 1 on every 
link. The dual action then simply reduces to that of a 
global Ising model for the Vi with the dual coupling Kd- 



APPENDIX D: DUALITY OF THE MODEL 
WITH COMBINED [/(I) AND Z2 INVARIANCES 

In this Appendix, we will perform a duality transfor- 
mation on the chargon-spinon action S = Sc + Ss + Sb 
derived in Section ^ to work instead with vortex variables 
instead of the chargons. For simplicity, we will restrict 
ourselves to situations with an integer number of elec- 
trons per unit cell. In this case, the Berry phase term 
Sb is independent of the chargon phase field (j)i . In Sec- 
tion VI, we will provide the generalization necessary to 

All 



handle non-integer number of electrons per unit cell, 
of our transformations will focus entirely on the term in 
the action involving the chargon variables. This is simply 
a chargon hopping term: 



Sc^-J2a,,{t,b*bj+c.c.), (Dl) 

^ 2ieC0s((/., -(/., + ^(1- ay)). (D2) 



<ij> 



Note that in the absence of aij^ this is just the action for 
the three dimensional XY model. The duality transfor- 
mation for the 3DXY model is standard - here we will 
generalize it to include the Z2 gauge field aij . 

Consider the partition function obtained by integrat- 
ing over the chargon fields in the above action: 



Zholi'^] = 



2Tr 



(D3) 
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As with the duahty transformation of the XY model, it 
will be convenient to work with the Villain form of the 
action 

S[(t), J,cr]= Y^ njfj/2 + iJ^j{(j)^ - cj)j + |(1 - CTy )), 



<y> 



(D4) 



where Jy are integer valued fields that live on the links of 
the lattice, and are to be summed over in the partition 
function. As usual, this is strictly justified only in the 
limit tc « 1 when tc — exp(— k/2), though we do not 
expect any modifications to the physics by relaxing this 
assumption. The Jij have the interpretation of being the 
total conserved electrical current on any link. This can 
be made more explicit by performing the integrals over 
(f)i which imposes the current conservation condition 



A • J = 0. 



(D5) 



The symbol on the left hand side is the lattice divergence 
of the link variable Jij. We proceed, as usual, by solving 
the current conservation condition by writing 



27r J,- 



Ax 



(D6) 



The quantity a lives on the links of the dual lattice, and 
is constrained to be 27r times an integer. The right hand 
side is the lattice curl of this variable a on the plaque- 
tte of the dual lattice pierced by the link < ij >. The 
chargon action now takes the form 

^[«' '^l = E Ft(^ X a)2 + ^ ^ (A X a)(l - a,,). 



^ 87r2 ' "' 'A . . 



(D7) 



Here the first term is a sum over plaquettes of the dual 
lattice, and the lattice curl in the second term is on the 
plaquette pierced by the link < ij >. Now note that 
as Uij = ±1, the exponential of the second term can be 
written 



n(-i)(^)(- 



It is useful now to separate the integer ^ into its even 
and odd part by writing 

a = 27r(2A + s), (D8) 

where A is an integer and s = 0, 1. Then, we have 

(D9) 



where the product inside the brackets denotes the prod- 
uct over the links of the plaquette of the dual lattice 
pierced by < ij >. We now define 



Ai»j = (-1)' = 1 - 2s. 



(DIO) 



Note that Hij lives on the links of the dual lattice and 
takes values ±1. The product above can then be written 



exp I 



5 (i-n.)(i --.)). 



Note that /x satisfies 



_i r»j 



(Dll) 



(D12) 



where the plaquette product on the left hand side is on 
the plaquette of the dual lattice penetrated by the link 
< ij >. Thus, the conserved Z2 charge current deter- 
mines the flux of ^. 

The action now becomes 



^ = E 



^ Ax 2A+^—^ 
87r2 I I 2 






Scs, (D13) 
(D14) 



At this stage, A is constrained to be integer-valued. We 
impose this integer constraint on A softly by adding a 
term 

-tyYcos{2nAij). (D15) 

Here the sum is over the links of the dual lattice. The ac- 
tion can now be rewritten in terms of a = 27r [2 A H — j^ j : 

S = Sy + Sa + Scs, (D16) 



u-i^ f^v / ^ 



[iij COS 



K 



^« = Ei^(^x«)'- 



Btt 



(D17) 
(D18) 



It is convenient to extract a "matter field" from the Oy 
by letting 

fly ^ fly +2(0,-0,). (D19) 

This changes Sy to 



Sv == -it, ^ fiijcos [di- 6j + -^ 



(D20) 



<y> 



but leaves all the other terms unchanged. The field e*^* 
may be interpreted as an -^ vortex creation operator. 
Several symmetries of the action above are apparent. It 
is invariant under a local U{\) gauge transformation 
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Or + A,; 



A,,; - A, 



(D21) 

(D22) 



This is standard in the dual vortex description of XY 
models in three dimensions. However the action has an 
additional Z2 gauge symmetry under which 



„»e. 



js. 






(D23) 
(D24) 



with ei = ±1. This Z2 gauge symmetry is actually dual 
to the one in the chargon-spinon action. Note that the 
action describes the vortices e'^' minimally coupled to the 
fluctuating U{1) gauge field a , and also to the fluctuat- 
ing Z2 gauge field fi. The field /i is in turn coupled to 
the field a by the term Scs- 

This completes the duality transformation to the vor- 
tex description. Adding together the spinon action and 
the Berry phase term Sb gives the full dual action of 
Section 0. 
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